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FOREWORD 


This is the Final Report on IIT Research Institute Froject K231, 
"Fundamental Study of Debris Hazards," conducted for the Defense Atomic 
Support Agency, Washington, D.C., under Contract DA. 49..146-XZ-097. 

All work done under the initial contract on the basic study of debris hazards 
which provided for supplemental research to analyze data and develop 
expressions to facilitate the formulation of debris damage predictive schemes 
for several specific situations done under Modification No. 1, are reported. 


Work performed on ''DANNY BOY TASK II - Debris Investigations, " 
under Contract Modification No. 2, was described in the preliminary and 
final test reports, ''Throwout Study of an Underground Nuclear Detonation, '' 
published by the Department of Defense - U. S. AEC as PORO1814 (ITR~1814) 
and POR-1814 (WT-1814). 


Work done on ''DANNY BOY TASK I - Pre-Shot Predictions," also 
under Contract Modification No. 2, was reported as a separate report en- 
titled, "Hydrodynamic Analysis for a Buried Underground Nuclear Explo~ 
sion."" Work performed on ''DANNY BOY TASK III - Post:-Test Crater 
Analysis,'' under Contract Modification No. 3 will also be submitted as a 
separate report. 


The ccoperaticn and assistance of the Armed Services Explosives 
Safety Board, (especially Mr. Russel G. Perkins, Chief of Explosives 
Branch) in allowing project engineers use of their explosions files 1s 
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IITRI personnel contributing to these studies include E. B. Ahlers, 
D. I. Feinstein, B. Gain, P. C. Hermann, J. Lukes and Dr. K. E. McKee. 
Mr. C. A. Miller conducted the analysis of horizontal motion of debris 
particles under the influence of the blast winds. Mr. R. L. Barnett 
developed the analysis of the motion of tree debris. 
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Project Engineer 
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ABSTRACT 


Predicting the mechanical debris, associated with nuclear detonation, 
stemming from several sources (blast-induced flight of structural fragments, 
pickup of material from the ground by blast winds, and crater throwout) is a 
significant problem. Designers of hardened sites are concerned with mater - 
ials which may accumulate atop silo doors or damage vulnerable above-grade 
antenna systems. Military operations are concerned with the hazards to 
troops and equipment from the use of nuclear demolition and tactical devices. 
Potential users of nuclear excavating devices are concerned with hazards to 
personnel, utilities and equipment. 


This rerort describes the collection and analysis of data on various 
aspects of debris formation and dispersion, and examples of data utilization 
in estimating debris environment in several situaticns. The approach in the 
study was to collect extensive data from past experimental and analytical 
investigations bearing on debris formation and dispersion. By further study 
the most meaningful formats were summarized to be used as inputs to 
approximate solutions for debris environment predictions. 


An extensive regression study of several hundred HE incidents, 
accidental and experimental, is made to relate the maximum range of debris 
to explosion parameters and crater dimensions. Results showing consis- 
tency with the limited available nuclear data relating to crater throwout are 
also presented to describe the nature of the debris distribution function in 
general terms. 


Fragmentation data from HE events and laboratory experiments are 
used to indicate the nature of fragment.-size distributions from structural 
demolition. 


An analytical study of the motion of debris fragments caused by 
blast winds considers debris trajectories for various times of structural 
failure, fragment sizes, positive and negative phase winds, and initial 
elevations of the fragment. 


Specific estimates are made of the debris hazards to troops of fly- 
ing tree limbs in the proximity of forest stands, the vulnerability of troop 
personnel to throwout debris from cratering and streambed charges, and 
the debris environment about hardened antenna systems. 


Useful estimates of debris environment can be made for many tar: 
geting situations with data contained inthis report. Refinement of data is 
certainly essential, especially in experimental definition of fragmentation 
patterns of ideal structural elements, and in definition of crater lip contours 
near their extremities, i.e., the throwout debris within and beyond the 
extremities of the lip. 
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CHAPTER ONE 


INTRODUCTION 


While most nuclear weapons effects -- radiation, thermal, air 
blast, ground shock, and radioactive fallout -- have been investigated ex- 
tensively since 1946, little study has been devoted to the creation and dis- 
tribution of structural debris. This stems in part from the fact that under 
earlier concepts of vulnerability, yields of weapons and the accuracy of 
their delivery, the design hardness levels under consideration were such 
that structural debris was not a prime hazard. Developments in the yield 
of weapons, delivery accuracy, and design hardness levels make the effect 
of flying debris a serious consideration. With the utilization of tactical 
nuclear weapons by field troops, knowledge of debris environment is an 
important consideration in the deployment of troops. Design specifications 
for survival of hardened retaliatory weapons sites shortly beyond the edge 
of the plastic zone of the crater require a knowledge of the deposition of 
crater throwout material. The siting of communications systems presents 
particularly serious problems since they are especially vulnerable where 


substantial flying debris arises. 


The investigation was initiated to conduct a series of analytical 
and experimental studies of the behavior of debris. The following tasks were 


included: 
Fundamental Studies of Debris Behavior 


Measures of the maximum fragment distance were 
developed from HE data and compared with nuclear 
results. The shape of debris distribution functions 
was studied in detail. Experimental and analytical 
work on fragment size-distribution was revised. 
An analytical model for the transport of debris 


under blast wind loading was developed. 
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@ Specific Vulnerability Studies 
Several vulnerability studies were included, both as 
applications of the collected data, and to provide 
measures of vulnerability to flying debris under 
targeting situations of considerable interest and 


significance. These studies included: 


® Vulnerability of field troops from branch wood 


within or near forest stands. 


e Vulnerability of field troops to crater-emanated 


debris from very-low yield cratering charges. 


e Vulnerability of field troops to crater-emanated 


debris from very-low yield stream bed charges. 


e Vulnerability of antenna systems crater-emanated 


debris from high-yield nuclear bursts. 


e Crater Throwout Study of an Underground Nuclear 


Detonation 


This task was conducted as DANNY BOY Project 1.5 
with findings reported under separate cover (Ref. 1, 2). 
DANNY BOY Project 1. 5 involved the following activities: 


e Compilation of data relating the initial and terminal 
positions of a series of more than 1100 "ideal" 
objects (steel plates, spheres, cylinders and cubes; 
wood cubes and boards; and common brick) emplaced 
on the ground surface and in drill holes in the crater 


zone prior to the shot. 


e Compilation of data on the distribution of natural 
throwout debris beyond the limit of the crater lip, 
i.e., beyond the ground range where the ground 


surface is completely obscured by debris. 


e Compilation of analysis and plotting of the data in 
various manners to describe the behavior of crater 


throwout for this deep underground burst. 
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ll Report Organization 


This report is organized into seven chapters. Each of the six 
chapters following the Introduction concerns a pertinent aspect of the over- 
all debris problem. Thus, the report is a series of related, but self- 
contained, studies rather than a continuous-reading document. The contents 


of the report are summarized by chapters as follows. 


Chapter One, Introduction 


The first chapter provides a general introduction to the debris 
problem. The major conclusions drawn from the investigation are presented. 


Finally, recommendations for additional research are made. 


Chapter Two, Debris Characteristics of High Explosive 


and Nuclear Detonations 


Previous nuclear weapon effects tests have produced little data 
concerning the formation and distribution of structural debris, with the 


following exceptions: 


Project 4. 5 of Operation JANGLE included measure- 
ments of the debris from airport-type runways and 
reinforced-concrete wall panels erected over the 
crater zone of an underground burst. 


Project 33. 2 of Operation PLUMBBOB studied the 
behavior of debris (including window glass, military 
debris, gravel, stones, and spheres) in response to 
air blast at various ground ranges. Findings of this 
investigation had not been published at the time of 
this investigation. 


Other Operation PLUMBBOB studies were concerned 
with fragments of biological interest, suchas glass 
splinters capable of penetrating abdominal walls. 


DANNY BOY Project 1.5, a study of crater 
throwout, was recently completed. 


SEDAN, a recent nuclear event, included an experi- 
mental investigation of crater throwout. 
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While debris information from past nuclear tests provided only 
limited data for analysis, extensive measurements were found to be avail- 
able in reports of planned and accidental HE detonations. These data were 
collected, studied, and plotted in various manners to define general patterns 
in debris behavior. Where possible, 'check points'' from nuclear experi- 
ence are introduced to indicate their consistency with HE results. Hundreds 
of explosion reports were reviewed in this task. Data from a series of 
more than 200 selected explosions were used to obtain expressions relating 
maximum debris distance to equivalent yield, and to an estimate of the 
explosive impulse. A smaller series of explosions was used to derive 
expressions relating maximum debris distance to crater dimensions. 
Maximum debris distance is correlated with explosion parameters and 
crater dimensions using the method of least squares; scaled and unscaled 
relationships are developed using both linear and quadratic relationships. 

In each regression line so obtained, the standard error and correlation 
coefficient have been computed as a measure of closeness of fit. Figure 1.1, 
which expresses the quadratic correlation between maximum debris distance 
and equivalent yield, and Fig. 1. 2, which shows the correlation between 
scaled maximum debris distance and scaled crater volume,are typical of 
these regression results. Note particularly, the consistency of JANGLE U 
results with the HE findings. The much lower position of the DANNY BOY 
results can be explained on the basis of the very deep burial of the device 

in this event, which resulted in trajectories with pronounced vertical 
components in the ejecta. Actually, it may be argued that the HE detona- 
tions are more like the buried nuclear burst than surface nuclear bursts, 


because of the absence of substantial blast winds. 


To define the probability of personnel or equipment being hit by 
missiles, and, for the shorter ground ranges, the quantity of material likely 
to be deposited atop a silo door, it is necessary to estimate the distribution - 
of fragments within the maximum ground range. Several approaches to this ase 
problem were pursued. A theoretical model for fragment distribution from 
ground zero to maximum debris distance was derived by assuming wall and 
roof panels to fragment into equisized fragments upon detonation of a near- 


by line charge. Comparison of this model with several explosions shows 
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that the shape of the curves is similar (Fig. 1.3). Next, debris distribution 
patterns of a series of six ordnance explosions were plotted in various mat- 
ters to note their similarity. Thirdly, debris distribution of an ordnance 
structure (involving over 30,000 recorded fragments with a total weight of 
about 43 tons) was plotted in detail. Contrary to expectation based on drag 
effccts, the large fragments from this explosion did not travel as far as the 
s:naller ones (Fig. 1.4). This is probably because fragments subject to 
forces Sufficient to cause large acceleration were also subject to a greater 


degree of breaxup. 
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Chapter Three, Fragmentation Experimental Observations 


Since energy levels of the flying fragments are a measure of 
their potential to penetrate shields, impart shock loading to equipment, or 
incapacitate personnel, it is well to define the expected size-distribution of 
fragments. No definitive experimental investigation of structural fragmenta- 
tion has yet been made, even on idealized structural elements, by which the 
fragment-size distribution is related to structural strength and loading 
parameters. For this reason an attempt was made to collect and summarize 
past experimental work on fragmentation of materials and structures, and 
to describe the general fragment-size distribution patterns that have been 
observed. Five such investigations are described; 

The British Coal Utilization Research Association 

studies on coal breakage from random forces, 


Safety in Mines Research Establishment (of Great 
Britain), research on explosively detonated stone 
blocks. - 


Stanford Research Institute model tests on frag- 
mentation of reactor containment structures 
from internal explosions, 


The Pantex Ordnance Plant detailed fragment counts 
from the planned explosion of a reinforced concrete 
ordnance structure, 


Project 4. 5 of Operation JANGLE studies of fragment 

size distributions from reinforced concrete wall 

panels erected over the crater zone. 
Experimentation has shown that the higher the loading on the source material, 
the smaller the fragments produced, and that a wide range of fragment sizes 
are produced by any loading. An "Ideal Law of Breakage" which shows 
excellent fit to the experimental data has evovled from coal-crushing 
investigations. The extensive data from over 30,000 concrete fragments 
in the Pantex Ordnance Plant event allowed a detailed study of the fragment- 
size distribution. It is interesting to note that only about 3 percent of the 
fragments recorded in this event (above l-ounce in size) weighed more than 
three pounds, but that these accounted for nearly 75 percent of the total 
weight of all fragments. This tends to support the hypothesis that for 


many problems involving impact of fragments, there may be an optimum 
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fragment-size for purposes of design predictions. Figures 1.5 and 1.6 


show these distributions. The only major structural fragmentation study 


conducted on nuclear detonations was the test of wall panels on the JANGLE U 


event. Size distributions of the larger fractions were plotted as part of 
that project and it was noted that the JANGLE data did not preclude the 
possibility that the concrete fragment-size distribution caused by the 


underground nuclear shot followed the same pattern as mined coal or ore 


in a crusher. 
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Chapter Four, Fragmentation, Analytical Considerations 


Prediction of the number and size into which a structural 
component fragments, when subjected to blast loading, is totally beyond 
the current state-of-the-art of fracture mechanics. Thus, a relatively 
simple mathematical model was formulated to be used to predict debris 
formation, which can be adapted to include advances in fracture mechanics. 
The model selected for this study, is a simple extension of the single-degree- 
of-freedom structural model long used to analyze the elastic-plastic 
response of structural components subjected to blast loading. The response 
of the single-degree-of-freedom system with an elastic, perfectly plastic 
spring is determined. In particular the mass, velocity, and the time when 
the mass reaches a certain displacement, called the fracture displacement, 
are of interest. The magnitude of the fracture displacement and the 
description of particles formed when this displacement is reached must 


await further analytical and experimental developments in fracture mechanics. 


Chapter Five, Debris Transport by Blast Winds 


The motion of a particle acted on by the nuclear blast winds is 
analyzed on the assumption that initial conditions of the particle motion are 
known. The model used assumes that the force acting on the particle is 
proportional to the square of the relative velocity between the particle and 
air. Blast parameters are assumed constant over the range of travel of 
the debris, and it is further assumed that the apparent lengthening of the 
positive phase duration due to the debris motion in the direction of shock 
propagation can be handled by a simple adjustment of positive phase duration. 
These assumptions are necessary to reduce the equation of motion for 
any debris particle to a one-parameter nonlinear differential equation. 
Without making these two assumptions, it would be necessary to treat each 
weapon yield and placement as a separate protlem and no general obser- 
vation could be made regarding debris behavior. The equation of motion 
is numerically integrated and results are obtained for a wide range of 
overpressures and particle sizes. It is possible to use the results directly 


to determine distance-time plots for any debris. The effect of negative 
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phase winds on particle motion is also studied. While a complete 
description of the negative phase wind is not available, reasonable 
estimates of this wind can drastically change the motion of the debris 


particle. 


Chapter Six, Vulnerability of Field Troops to Tree Debris 


The study of the vulnerability of engineer and field troops to 
hazards from tree debris is one of three specific debris problems of in- 
terest to the Office of the Chief of Engineers, U.S. Army. The objective 
of this study was to define a ''safe distance" for positioning troops to avoid 
casualties from falling trees and limbs. The safe distance may actually 
fall inside or outside the forest and both cases are represented. Since 
any individual tree may fail at the base under blast loading: tree height 
is considered the minimum Safe distance in any situation. Making certain 
simplifying assumptions (i.e., zero strength tree limts, plane blast wave 
loading, unobstructed trajectories, and that personnel struck by tree limbs 
are certain casualties),tree limbs are followed in their trajectories from 
the time of shock arrival to their impact with the ground. Results show 
that for the lower yields (1 KT, for example) trajectories become vertical 
early. A uniform translation of all branches is thus obtained, the area in 


front of the forest up to the ''safe distance" being similar in appearance 


to that of the forest floor after all branchwood was allowed to drop vertically. 


For the higher yield weapons (20 MT, for example) trajectories terminate 
before they become vertical. Results are similar with the exception that 
the highest branches of the first few rows of trees pile up in a lower 
density than those following closer-in trajectories. Results of this study 
are summarized in the following table. which lists safe distances in terms 


of tree height, weapon yield, and overpressure levels. 


Chapter Seven, Vulnerability of Field Troops to Throw Out Debris 
from Cratering and Stream Bed Charges 
Methods for estimating safe distances for positioning troops in 


the proximity of very-Jow-yield cratering and stream-bed charges, based 


on debris criteria; were studied. As with the preceding tree debris vulner- 
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SAFE DISTANCES TO PREVENT CASUALTIES FROM TREE DEBRIS 


Safe Distance “Falls a 
Inside of Forest Safe Distance from Ground Zero Safe Distance from Forest 
for Various Overpressures for Various Overpressures 
at the Front of the Forest, at the Front of the Forest, 
from Ground Zero _{yd) _(yd) 


Seta pel Bees | ipa ee) eee tel TO pai] 15 pat [20 pai 
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ability problem, these two cases were studied to supply data needed by the 
Office of the Chief of Engineers, U. S. Army. Numerous reports on debris 
from HE cratering tests were studied to obtain measures of the relationship 
between debris distance, weapon yield, depth of burst, fragment size, and 
soil characteristics. The U. S. Geological Survey report on a series of 
HE cratering tests in basalt (Area 18 at Nevada Test Site) provided data 

on relationships between the first four of these factors (Ref. 24). Using 
these data, nomographic procedures were presented for estimating average 
debris distance for various fragment sizes in basalt, for various combin- 
ations of yield and depth of burst. An asymmetry factor is introduced to 
obtain maximum rather than average debris distances, based on measure- 
ments of the ray-like patterns in the USGS study. Using data from the 
Panama Canal series of cratering tests in various media (Ref. 29), distance 
ratios are plotted obtaining estimates of debris distances in media other 
than basalt, and for the stream-bed charge. Debris distances found by 

this method, using DANNY BOY explosion parameters, are consistent with 


photographic observations of the debris deposit from this event. 


1.2 Summary and Conclusions 


Various aspects of the debris problem -- fragmentation charac- 
teristics of materials and structures, transport of the debris particles by 
the blast winds, and the ultimate distribution of material resulting from 
these factors -- are studied in this investigation. Both analytical and 
experimental studies of the phenomena are considered. In some cases, 
behavior is defined empirically from collected historical data, in other cases, 


by analytically derived expressions and exhibits. 


Fundamental data are collected to provide a basis for initial 
estimates of the nature of the debris problem in specific targeting situations. 
Historical data (HE detonations) are used to derive empirical expressions 
correlating maximum debris distance with equivalent TNT yields, impulse, 
and crater dimensions. These serve to make initial approximations 
of the limiting ground range of debris problems. Data compiled on relative 
distribution of debris particles can be used to make initial estimates of 


fragment quantities at various ground ranges. Estimates made from these 


SECRET 


15 


een nn anne 


SECRET 


data are, at best, crude. For the most part they are based on HE data 

and nuclear - HE equivalence. Use of the latter is questionable -- expecially 
since extreme differences in blast winds exist between nuclear and HE. As 
with most HE results, the experimental findings include great degrees of 


scatter. 


The existence of definite debris patterns is graphically demon- 
strated in this report. The only pehomenon lending itself to analytical 
treatment is the transport of debris particles by the blast winds. A 
relatively complete treatment of the phenomenon is made. Fragmentation of 
materials and structures is reviewed in detail -- but this study simply 


cannot be carried very far without considerable experimentation. 
The following conclusions are drawn from this study: 


(1) Data in this report can be used to determine debris 
distribution under a variety of targeting situations. 

(2) In many situations initial estimates of the severity 
of the debris problem can be made from the data 
compiled. 

(3) Progress has been made in providing data and 
methods for estimating the existence and severity 
of debris problems, but the approximations are crude. 

(4) The various debris phenomena (fragment-size distri- 
butions for materials and structures, debris disper- 
sion patterns, and variations with explosion parameters) 
follow characteristic patterns which can be established 


experimentally. 


(5) The feasibility of predicting debris effects of nuclear 
explosions has been demonstrated by examples. Be- 
cause random and/or uncontrollable (and unpredictable) 
factors influence debris behavior, prediction of these 


effects will never have a high degree of precision. 


(6) Debris prediction methods can be greatly improved by 
experimental data which can be developed on future full- 


scale nuclear test programs. 
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Problems exist which cannot be solved without further investi- 
gation, primarily experimental. Knowledge of fragmentation of actual 
structures, or even structural elements, is sub-minimal. Little is known 
of fragment quantities and the fragment-size distribution which can be 
expected from structural elements, and even less is known of the variations 
in the characteristics which would be caused by different levels of impulsive 
loading. Likewise, the time of failure of structural elements is a completely 
unknown factor. The analytical treatment of fragment transport by blast 
winds, included in this report, requires experimental verification. The 
problem of crater throwout debris and the resultant buildup of the crater 
lip is significant in terms of the mass of material likely to be accumulated 
atop hardened sites. The work on crater throwout debris begun on 
Project DANNY BOY under this contract should also be expended to cover 


variations in parameters. 


1.3 Recommendations 


On the basis of the findings of work performed under this 
contract, it is recommended that further analytical and experimental 


investigation te made in the following areas. 


(1) Further study of the data accumulated in this report, 


aimed toward codifying estimating procedures. 


(2) Further analysis of the accumulated data from the 
crater throwout debris study from the DANNY BOY 


event. 


(3) Inclusion of crater throwout debris studies in possible 
forthcoming tests, to study effects of variations in 
parameters -- other yields, depths of burial, and 


soils, 


(4) Experimental investigation of fragmentation of photo- 
type wall and roof panels -- studying time of failure, 
fragment-size distribution, fragment dispersion, and 
effects of variations in the loading pulse. This work 
could be performed under a combination of HE and 


full-scale nuclear testing. 
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CHAPTER TWO 


DEBRIS CHARACTERISTICS OF HIGH EXPLOSIVE 
AND NUCLEAR DETONATIONS 


Little data on formation and dispersion of structural debris have 
been collected from full-scale nuclear tests. Since debris hazards were not 
initially considered a major problem, tests of structures and structural 
elements did not include measurements of fragmentation or debris transport. 
One exception was Project 4.5 of Operation JANGLE, in which reinforced 
concrete wall and runway panels were erected above the expected crater zone 
and debris measurements taken (Ref. 5). The vulnerability of parked air- 
craft was emphasized in this study, thus major interest was centered about 
the maximum ground range of aircraft damage from blast and from debris. 
A second exception was Project 33. 2 of Operation PLUMBBOB in which the 
behavior of missiles (including window glass, military debris, gravel, 
spheres and native stone) emplaced at various ground ranges was studied. 
Findings of Project 33. 2 Operation PLUMBBOB had not been published at 
the time of this investigation. Other debris studies on Operation PLUMBBOB 
involved such fragments as small glass fragments capable of penetrating 
abdominal walls, which are of biological interest but which cannot be re- 
garded as structural debris. More recently, crater throwout debris studies 
included in underground nuclear test programs on DANNY BOY and SEDAN 
have investigated objects emplaced on the ground surface and within the ex- 
pected crater, and their post-shot locations relative to original positions. 
These crater throwout studies involved deep-buried shots in which the 
debris was not influenced substantially by blast winds, as it would be witha 


surface or shallow-buried shot. 


From the beginning of this investigation, it was apparent that frag- 
mentation and debris dispersion data from past full-scale nuclear tests would 
provide only limited information for analysis. It became obvious, however, 
that an extensive body of data was available in explosion reports from HE 
events -- both planned and accidental -- and that these could be studied to 


define general debris behavior patterns. The limited nuclear data could 
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then be introduced as ''check points" on the HE studies. Several hundred 
selected reports on HE detonations were studied to determine the 
relationship between the maximum debris distance and various explosion 
parameters; these were then used to define the expected limit of the debris 
hazard. The least-squares method was used to correlate maximum debris 
distance with equivalent yield, computed values of impulse, and crater 
dimensions. It is of particular interest to note that comparisons of maxi- 
mum debris distance from the JANGLE and DANNY BOY studies are consist- 
ent with the HE results. 


In addition to defining the maximum debris distance, it was 
necessary to estimate the distribution of debris fragments at distances 
shorter than maximum. Quantitative data are needed for any estimate of 
probability of equipment or personnel being hit by flying fragments in the 
debris zone. Several approaches to this problem were taken. A theoretical 
model for the distribution of fragments from ground zero to maximum debris 
distance was developed by assuming structural wall and roof panels to fail 
into equi-sized fragments, with radial initial velocity vectors, upon detona- 
tion of a line-charge placed near them. Comparison with debris-distribution 
patterns of several explosions show that curves are of similar shapes. Next, 
the debris-distribution patterns of a series of six HE detonations involving 
ordnance structures were plotted together in various manners to note their 
similarity. Thirdly, debris distribution from a reinforced concrete ordnance 
structure was plotted in various ways. The availability of data, from this 
event, on the weights and terminal locations of over 30,000 fragments with a 
total weight of about 43 tons permitted highly detailed treatment of debris 
distribution. An interesting observation of this study is that the large frag- 
ments did not travel as far as the smaller fragments. This seemingly 
contradicts expectations based on air drag effects, but the comparison is 
inappropriate since the particles subject to forces sufficient to cause large 


motions are also subject to a greater degree of fragmentation. 


2.1 Application of HE Data from Planned 


and Accidental Explosions 


Past HE detonations, both planned and accidental, are an exten- 
sive source of data on the behavior of structural debris. Explosion reports 


customarily include some or all of the following information. 
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Total weight of explosive involved 

Weight of explosive exploding at one time 

Kind of explosive 

Source structure 

Crater dimensions 

Maximum debris distance 

Distances at which individual fragments were found 

Distances at which various degrees of structural 

damage occurred 
Several detailed compilations of explosion reports have been made (Ref. 6, 7, 8}. 
Explosion reports make it possible to plot debris distances -- frequently 
only the maximums -- against various explosion parameters. Statistical 


measures of the consistency of these data can also be made. 


Although the validity of debris data from individual explosion 
reports is limited in certain respects, the availability of many reports does 
permit averaging the data and constructing graphical and statistical mea- 
sures of debris behavior. Plots of maximum debris distance in terms of 
equivalent explosive weight which permit estimates of expected debris 
ranges are included here. While scaling from HE to nuclear is generally a 
questionable procedure, these plots can be, and in fact have been, applied 
in determining whether or not debris problems are likely to exist at various 


ground ranges. 


Certain limitations in the use of data taken mostly from reports 
on accidental explosions are apparent. Extreme variations exist in the range 
of explosion environments. Explosion source structures include steamers, 
freight cars, munitions plants, warehouses and igloos. Some of the explo- 
sions occurred in the open field. The relative degree of confinement afforded 
by these structures certainly affects the loading on the flying fragments. 

The amount of material available for fragmentation is a function of the source 
structure and, for larger explosions, of the surrounding neighborhood as 
well. In some cases, explosion data may be questioned as to the accuracy 

of the weight of explosives involved in accidental explosions, whether high 

or low order detonation took place, and the quanity of explosives involved 


in individual detonations of a series of consecutive explosions. Also, the 
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trajectories of fragments from HE detonations will range from the horizontal 
to the vertical, the farthest-flung may have initial trajectory angles of around 
45° from the horizontal. This behavior is unlike fragment trajectories re- 
sulting from a structure failing under loading from the plane blast wave of a 


nuclear explosion. 
2.2 Measures of Maximum Debris Distance 


A series of 206 accidental and planned explosions was selected 
from those tabulated in references 6 and 8 and subjected to a regression 
study to develop estimating expressions for maximum debris distance -- a 
measure of maximum range of vulnerability to debris hazards. Selected 
explosions are listed in Appendix A. which tabulates maximum debris 
distances and the major explosion parameters. Maximum debris distance 
was correlated with equivalent weight of explosives and impulse; various 
scaled relationships between these functions were also studied. A total 
of 40 linear and 40 quadratic regression lines was determined by computer 
methods, together with standard errors and correlation coefficients for each 
line. (Results of this correlation study are included in Table 2.1 and Fig. 
2.3 through 2. 10.) 


2.2.1 Details of Regression Analysis 


The plots of maximum debris distanc: und explosion parameters, 
(Fig. 2.3 through 2.10) show a great degree of scatter among data points. It 
is quite apparent that this dispersion is, in all cases,far too great to permit 
visual location of any average line through the data points. Statistical devices 
must be used to describe the plotted relationships satisfactorily. Simple 
(linear) and quadratic correlations were therefore made, using least-squares 


methods. 


It should be noted that the computation of a regression line does 
not necessarily assure a functional relationship between the factors corre- 
lated; in this case it merely describes what has historically occurred. Thus, 
the computation and drawing of a regression line describing the correlation 
between maximum debris distance and equivalent weight of explosives does 
not inean that the two factors are functionally related by some physical law 


to the exclusion of such other considerations as configuration and strength 
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of structures and impulse distance. If, however, the association between 
maximum debris distance and some explosion parameter were found to be 
sufficiently close, it could be possible to estimate, with some calculable 

degree of accuracy, the maximum debris distance for other explosions on 


the basis of the observed relationship found in this sample of 206 explosions. 


The computed regression lines were derived by the method of 
least squares -- one of the most widely used methods of curve fitting -- 
which yields the equation of the line from which the summation of the 
squares of the deviations of all the plotted points is a minimum. Application 


of the least-squares linear regression line is demonstrated in Fig. 2.1. 


Derivations of equations for least-squares regression lines can 
be found in most engineering statistics texts. The linear least squares line 


of the form 

er Spat ya ho 
is obtained by solving the simultaneous equations 
+b), 2 X,= 2X 


2 12 2 1 


ro 
D> = 
aj2> X,+b),2% X, zX)X, 


Na, 


to obtain the coefficients al? and b> 


Similarly, the quadratic least-squares regression line of the form 
2 


"+b! Xj) +c' X 


i913 2 2 


is obtained by solving the set of three simultaneous equations: 
2 


t 1 ' _ 
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The advantages of deriving least-squares lines in this study are 
that (1) reproducible associative relationships are developed in accordance 
with accepted statistical procedures, and (2) measures of probable error 


and closeness of fit can be computed. The standard error is a measure of 
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X1.2 Unexplained Deviations 


Figure 2.1 Graphical Representation of Least-Squares 
Linear Regression Line 


X,, Xz = a set of bivariate data 


X, = arithmetic mean of X, 
Xi2 = least-squares regression line values of X, 
X)2 = deviations explained by regression of X, on X2 
Xi.2 = KX, - X,zg = residual or unexplained deviations 
N = number of items in sample. 
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the distribution of data points about the regression line, and is defined as 


(unexplained deviationel’ 


n=l ; 


Degrees of F reedom 


n 
2 
ou (X) - X,,) 


1 (Linear least-squares line) 


-_ 
Nw 
1 


(2. 2) 


(Quadratic least-squares line) 


The standard error about the least-squares regression line may 
be interpreted ina manner analogous to the standard deviation of a frequency 
distribution as represented by the normal curve. It is a measure of the 
distribution of data points about the least-squares regression line, and its 
interpretation is demonstrated graphically in Fig. 2.2. The upper chart in 
Fig. 2.2 shows the relative frequency of events producing data points within 
various multiples of the standard error about the regression line if the dis- 
tribution is normal. About two-thirds of all values would be within 
X,5+ 8) 2? and all but 27 out 10, 000 would fall within X,, + 3s 


12—- 12—""1,2 
lower chart is a more appropriate expression of this distribution for 


The 


vulnerability analysis. This chart shows the frequency of events occurring 
below the various multiples of standard error about the regression line. 


Its use is best shown by example. 


If we let D) equal maximum debris distance and W>: and W3 
the various equivalent yields, then the probability of a piece 

of debris being thrown as far as D, or beyond is 0. 14 percent 
for Wy) 2.9 percent for W> and 15. 9 percent for W Similar- 
ly, for a weapon of yield Ww, the probability of a piece of debris 
being thrown a distance equal to or beyond D, is 15.9 percent; 


and is 2.9 percent for D, and 0. 14 percent for D3. 
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Figure 2.2 Normal Distribution of Events 


about Least-Squares Regression Line 
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Standard error is expressed in units of the dependent variable 
and, as such, is not amenable to comparison with unlike quantities -- or 
moreover, with any general standard of closeness-of-fit. The correlation 
coefficient is a measure of scatter in dimensionless terms, and is defined 
as the square root of the ratio of the explained variation xt, to the total 


variation Ext, 
2 
2x12 
r= ay Aamo 
2x) 


The correlation coefficient can be computed for the linear case by the ex- 


pression 


2 
aS eae. meats 


The sign of r is positive for a regression line of positive slope, and nega- 


tive for a regression line with a negative slope. 


In general a correlation coefficient of +0. 90 or greater indicates 
high positive correlation, and between zero and +0. 10 indicates a low posi- 
tive correlation. Positive correlation coefficients below +0. 90 generally 
do not engender high degrees of confidence. A correlation coefficient of 


+0. 50 or thereabouts is decidedly marginal. 


2.2.52 Regression Study of Maximum Debris Distance 


The quantities of explosives involved in the 206 explosions studied 
ranged from 8 lb to 9, 000, 000 lb, providing more than seven cycles of data. 
The total weight of explosive materials involved was about 50, 000, 000 lb. 
Since the largest explosion in this compilation involved 9, 000, 000 lb (4.5 
kilotons of ammonium nitrate), extension of the regression lines to the 


magnitudes of nuclear weapon yields involves an extrapolation of one or two 
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orders of magnitude beyond the limits of plotted data. Their use in estimating 
maximum debris distance in the range of high-yield nuclear weapons requires 
extrapolation of three or four orders of magnitude beyond the limits of the 
plotted data. Thus, while there is no specific interest in the very small HE 
detonations, they are included in the total sample to extend the range of 

basic data -- while extrapolating several orders of magnitude beyond the 
limits of the data can be considered questionable under many circumstances, 


it is better to do so with seven cycles of data input than with a lesser amount. 


A total of 40 linear and 40 quadratic regression lines, in logarith- 
mic terms, was derived by computer methods. Each of the debris distance 


parameters, 


Log) 9 maximum debris distance (ft) 

3 
ed 
2/3 


Log), cube-root maximum debris distance (it! 


Log) 9 2/3-power maximum debris distance (ft 


) 


Log) g w!3_ scaled maximum debris distance ae 
| tonsay7 


2/3 


{ \ 
Log) 9 W -scaled maximum debris distance | at | 


tons c 
TNT 
was correlated against each of the explosion parameters: 


Logg equivalent yield (tons.ny7) 


Log 10 impulse 


lb-msec 
sq in. 


Log) 9 cube-root equivalent yield (tons.nnpp) 


, | lb-msec 1/3 
Logjg cube-root impulse ( —eqann” ) 


Logio 2/3-powerequivalent yield (tons. aA ) 


lb-msec 
sq in. 


2/3 
Po 


Log 10 2/3-power impulse ( 


lb~msec 
1/3 
TNT 


Log wl/3 scaled impulse 


in - tons 
Tou. W-'° seated deipulee: 2 meee 


lin - tonsany 
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The selection of associative relationships to be correlated was 
arbitrary. No analytical basis for a fixed relationship existed, and no 
justification for scaled relationships was apparent. By computer methods, 
the computation of a large number of correlations required little additional 
effort beyond that for the one or two obvious correlations, and provided 
the opportunity to check the relative closeness of fit for various approaches. 
The computer program outlined in Appendix B was written to type out 
coefficients for the regression lines, the standard error and the correlation 
coefficient. Regression study results are tabulated in Table 2.1. As this 
table shows, most correlation coefficients were between 0.51 and 0. 69, 
which is not a high positive correlation. No single correlation appeared 
significantly better than all others. The quadratic correlations were only 
marginally better than the linear correlations. The improved correlation 
coefficients for the quadratic cases can sometimes be regarded as suspect, 
in fact, especially if they yield expressions indicating a continually increasing 


slope with increasing explosion size. 


Several correlations were selected for plotting, together with their 


respective data points. These are presented in the following figures. 


Fig. 2.3 Maximum Debris Distance Vs Equivalent Yield 
(linear log-log scale) 


Fig. 2.4 Maximum Debris Distance Vs Equivalent Yield 
(quadratic log-log scale) 


Fig. 2.5 w)/3 scaled Maximum Debris Distance Vs 
Equivalent Yield 


Fig. 2.6 w!/3 scaled Maximum Debris Distance Vs 
Equivalent Yield 


Fig. 2.7? Maximum Debris Distance Vs Impulse 
(linear log-log scale) 


Fig. 2.8 Maximum Debris Distance Vs Impulse 
(quadratic log-log scale) 


Fig. 2.9 w/3 sealed Maximum Debris Distance Vs 
w!/3 scaled impulse 
(linear log-log scale) 


Fig. 2.10 w!/3 scaled Maximum Debris Distance Vs 
w!/3_Scaled Impulse 
(quadratic log-log scale) 


Text follows on page 51 
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REGRESSION LINES FOR MAXIMUM DEBRIS D *€ 


Type Logarithmic Form 
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Correlation Regression Line 


3 
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Table 2.1. | 
S FOR MAXIMUM DEBRIS DISTANCE IN TERMS OF EXPLOSION PARAMETERS 
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Maximum debris distance data from the JANGLE U shot and the 
DANNY BOY event are included in Fig. 2.3 to 2.6 to show their consistency 
with the HE data. Maximum debris distance for the DANNY BOY event was 
taken et 750 ft, the distance to a lobe of detached dust observed in post-shot 
aerial photographic observation (Ref. 9). The DANNY BOY event involved 
a 0.430-KT device buried at 110 ft. Maximum debris distance for JANGLE U 
was taken at 550 ft -- the farthest-thrown recorded fragment of the rein- 
forced concrete runways and wall panels built within the crater zone. The 
JANGLE U shot involved a 1. 2-KT device buried at 17 ft. These points are 
consistent with the HE results, especially when considered in relation to 
the regression line and the one-standard-error limit. The DANNY BOY 
event data compare less favorably with the HE regression lines than do the 
JANGLE U findings. The much greater depth of burial in DANNY BOY, 
where the device was placed at optimum" depth, accounts for this. The 
resultant trajectories observed in DANNY BOY had pronounced vertical 
components. The favorable comparison of the underground nuclear events 
with the HE results is not surprising. Because of the absence of extreme 
blast winds, the HE detonations may be more closely akin te the underground 
nuclear burst than to surface or above-grade nuclear bursts, as far as 


debris behavior is concerned. 


Lack of similarity between contained HE explosions and nuclear 
explosions under various siting and target conditions is obvious. The appli- 
cation of these curves to specific targeting situations can certainly be ques- 
tioned. The dearth of actual debris measurements from full-scale nuclear 
tests :s also realized. In view of these factors, the necessity of using avail- 
able data suchas this, tempered by engineering judgment, is apparent. 

Use of these charts will permit the making of order-of-magnitude estimates 
of maximum debris distance. Such estimates can at least indicate whether 
debris is a problem under various targeting conditions. A first-order es- 
timate of maximum debris range measured from ground zero can be obtained 
from the correlations of maximum debris distance and equivalent yield. The 
correlations of maximum debris distance with impulse can be applied to the 
problem of an individual structure located intermediate between ground zero 


and the target whose vulnerability is questioned, Consideration of the 
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standard error provides an indication of the likelihood of the target being 


within maximum debris distance. 


These charts, in themselves, do not provide any means for 
estimating the probability of the target being hit by debris, or estimating 
the number of pieces of debris likely to hit the target. They estimate 
maximum debris distance only-- the location of the furthest-thrown fragment. 
The problem of areal density of debris (the number of fragments per unit area) 


will be considered separately, 


Considerable difficulty was experienced in keeping within the 
accurate capacity of the electronic computer in these regression studies, 
The matrix solution to the simultaneous equations for the least-squares 
lines involves taking small differences between very large numbers which 
are summations of products or summations of powers of numbers. As these 
summations increase in size, the process taxes the accurate capacity of 
the computer. Although where necessary, recourse was made to double 
precision (16 digits) on the IBM 7090/401 computer, it was still not possible 
to ascribe a high degree of mathematical precision to the constants in the 
regression results. Thus we see that mathematically the correlation co- 
efficient for the quadratic case should not be lower than that in the linear 
case, or else the quadratic expression should have a zero coefficient to the 
squared independent variable. Regardless of this, we have observed that 
computed values of correlation coefficients were sometimes lower for the 
quadratic case, and that computed values of standard errors were frequently 
higher for the quadratic case. In these instances the differences in values 
of correlation coefficients and standard errors for the linear and quadratic 
cases were small -- actually marginal. The first digit of the standard error 
and correlation coefficient is good, the second digit probably merits some 
confidence, The constants in the regression lines are of unknown precision. 
The relative position of the lines among the scattered data points and the 
ratio of numbers of data points within and without the plus-and-minus-one 


standard error bands are reasonable. 
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2.2.3 Correlation of Maximum Debris Distance 


with Crater Dimensions 


Maximum debris distance was correlated with crater dimensions 
using data from the series of thirty-six explosions tabulated in Appendix C 
(Ref. 6). Results of the regression study are tabulated in Table 2.2, Linear 
and quadratic plots of maximum missile distance in terms of crater dimen- 
sions are presented in Fig. 2.11 through 2.16, while Wl oo getted relation- 
ships are presented in Fig. 2.17 through 2.22. JANGLE U AND DANNY BOY 
results are also plotted on these figures. Again the considerable depth 
of burial in the DANNY BOY event accounts for the position of these points well 


below the lower one-standard-error band. 
2.3 Debris Dispersion Patterns 


Theoretical models for distribution of debris about line charges 
are developed and compared with actual results from semi-contained explosions. 
Measured debris dispersion from HE detonations is studied. Comparisons 
of dispersion patterns for a series of six HE events of different magnitudes 
are made. Dispersion patterns for a single HE event are also studied in 


considerable detail. 


2.3.1 Theoretical Models 


Debris Dispersion for Line-Source Explosions 


Theoretical models of the dispersion of structural fragments in 
terms of the maximum missile distance were derived for flat roots, walls, 
and arches. These models were prepared for the purpose of comparing 
theoretical distributions with those from experimental and accidental explo- 


sions. 


Models developed to explain the dispersion of structural fragments, 


in terms of terminal ground range, are based on the following assumptions: 


Complete fragmentation of the structure 
into equal-sized fragments occurs 


Fragments follow ballistic trajectories, and 
air resistance can be neglected 
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All fragments have equal initial velocities, 

though other assumptions as to the initial 

velocity field may be made, 
These assumptions are obviously not true, but refined expressions either 
cannot be obtained for the general case or would unduly complicate the 
model, These assumptions are justified in development of a model for 
order-of-magnitude estimates of the relative densities of debris at various 


ground ranges. 


Figure 2. 23 represents a wall subjected to an explosive impulse 


from a line-source charge. 


Line Source of 
Charge 


Figure 2.23 Wall Panel under Explosive Impulse from Line 
Source Charge 
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In this figure, 


7 Tr 
1 cos @ 

ho, IR =: “ead 

cos @ ~ 

cos @6 
' 

Se = —+—_ inwhich the limit on @< tan”) 4, 

cos. @ 


and the fragment distribution in the wall becomes: 


Number of fragments 1 h' 


1 ; nes - 
< paca 
nit elevation angle be cos” @ with limit @< tan r 


Similarly, the geometry of the roof is shown in Fig. 2. 24. 


Line 
Source sy 


Charge 


Figure 2.24 Roof Panel Fragmented by Line-Source Charge 
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Here, 
r. dé 
dx = 2 
cos O 
h' 
2° Sin} (2. 4) 
t 1 
ee ee in which the limit on @2 tan"! = 


and the fragment distribution in the roof becomes 
i ' 
Te raements ot == See os with limit @> tan eeu : 
Unit elevation angle sin @ cos 6 r 


For the circular arch, shown in Fig. 2. 25, 


© Line Source Charge 


Fig. 2.25 Structural Arch Subject to Impulse Load From 


Line-Source Charge 


the distribution of fragment quantities becomes a constant. 


Number of fragments _ dx _ ie sae aa 
Unit elevation angle d@ 3 . 


The distance traveled by the individual fragments is represented in Fig. 2. 26, 


2 ines a eset” ett Sit 


Figure 2.26 Trajectory of Individual Fragment 
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in which, 
v= v_sin@ 
Vv fe) 
; 1 2 
h = Me sin @t «> gt 
and 
1 2 
h = Tr gt 
Therefore, 
3 sin @t = gt 
VG sin ® 
t = 
& 
Also, 
Yn = Yo £98 e 
Therefore 
ds 


v_cos 6t 
fo) 


se sin 8 cos @ 
ro) 
g 


and the trajectory distance parameter becomes 


doc sin 6 cos @. 


The fragment quantity parameters, — and ——— » are plotted 
cos @ 


against the fragment distance parameter sin @ cos @ in Fig. 2.27 for 

values of @ from 0° to 90°. The relative quantity of fragments landing 
within any circumferential sector can be found by integrating these functions 
over the appropriate limits of the function sin @ cos @. 

This integration is carried out as follows for walls: 

In Fig. 2.27 let 


1 
y = =a 
cos 8 
x = sin @cos @ 
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Maximuin fragment distances 
for walls, roofs, and cylinders 
are not necessarily equal, but 
would be a function or impulse 
distance to these structural 
elements. Hence the curve in 
the form presented here are 
not additive. 


Total Double- Values 
Curve for Walls 


0.6 


Fragment Distance Parameter (sin 9 cos 9) 


Figure 2.27 Theoretical Fragment Distribution Functions 
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Then 


_ dx 
dx = qo 3° 


(esin@-6:4 eos? 6)-G0 


(1-2 sin” 6) de. 


The relative quantity of fragments in strips centered about the line charge, 


(or circumferential strips centered about ground zero), then becomes: 


Q, =f yax 


@ 

=| (1-2 sin” 9} F 

, 9 
26 


1 cos 
a 
= tan ® - 2 (tan @ - 0) 
e 
1 
%2 oat 
= 20 - tan 0 with limit on @Stan"” — . 
@ r 
1 
Similarly, for roofs, let 
Sct 1 
2 sin 6 cos @ 
x = sin @ cos @ 
dx= (1-2 eine 8) de, 
“dO 
thus 
On = f vax 
ahs _(i-2 sin” @) 2 nine 8) 
dé 
8, “sin 9 cos © 

- 
= fos, tan @- 2 (-log, cos a] 

e 

1 
°2 -1 nh! 
= log. tan @ + 2 log. cos @ in which limit 6>tan ~ mas 
8 
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For arches: 


y = constant= k 


x = sin @cos 6 


ee (53 ata’ Oya; 


thus integration for arches yields a fragment quantity function as follows: 


e 


1 


Qn =f (1-2 sin®'6) do 


8 sin 20 
E - 2k is Ser saa 


Le ein-ze-|%2 
iy ae aes 


1 


Relative fragment quantities for circumferential bands about the 
point of burst, (assumed here to be at zero elevation at the center of the 
containing structure}, are obtained by equating Q over the appropriate 
values of %) and 8,- For wall panels 8, is 0° and 8, the angle included 
between surface zero and the plane through the line charge and the wall-roof 
intersection. For actual structures Q must be corrected to account for 
the different material quantities in wall and roof panels, and for differentials 
between normal distances between the line charge and the wall and roof 
panels. Thus, for the structure pictured in Fig. 2.28, the relative number 


of fragments falling within any band By - R, about the structure would be 


J Containing Structure 


Assumed 
Location of Line 
Charge 


Figure 2.28 Typical Structure Subjected to Fragmentation 


SECRET 


by Internal Detonation 


75 


SECRET 


Q- = t, rQ  +tphQ, 


ar 90° 
= t r 120 - tan @ +t,h |log tan@+2log cos @ 
w ° R € € 
0 @ 
A 
where 
ty = wall thickness 
t. = roof panel thickness 


Comparison of Ideal and Actual Debris Distribution 
from Semi-Contained Explosions 

A comparison of the ideal fragment dispersion functions with 
three actual HE explosions is made in Fig. 2.29. In making this comparison 
the quantities tie 8: ts and h were dropped because of lack of data fully 
describing the containing structures. Furthermore, the functions Q. and 
Q,, in Fig. 2. 29 are plotted for all values from 0° to 90°, since the actual 
value of aN was not known. Actually, this means that Q. and Q, are 
plotted for panels of infinite length. The theoretical curves Q. and Q. 
were plotted by equating (sin @ cos 8) Aa = 1.0, computing Q.. and Q, 
for 10 percent increments of (sin @ cos ) ee and plotting values of Q. 
and Qy at the midpoints of these sectors. Curves for the actual explosions 
were established by equating the maximum debris distance to 1.0, and 
plotting actual reported fragment counts for 10 percent increments of the 
maximum reported distance at the midpoints of the circumferential bands. 


The following explosions are included in this analysis. 


Wolf Creek 1944: (Ref. 10) 


Accidental Explosion. 

4, 800 lb of ammonium nitrate exploded in a building. 

305 fragments were recorded within a maximum 
range of one mile. 


Kankakee, 1943: (Ref. 10) 


Accidental Explosion. 

5, 300 lb of Bi-Oil exploded in a building. 

160 fragments were recorded within a maximum 
range of 3,750 ft. 
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Pantex Ordnance Plant, 1960: (Ref. 11) 


Experimental explosion. 

2,000 lb of HE in encased warheads detonated 
inside U-shaped bay of reinforced concrete 
ordnance structure having one-foot thick walls. 

About 31,000 concrete fragments, with a total 
weight of about 85,000 lb, were recorded within 
a maximum range of about 1,450 feet. 


The plotting of data from these explosions assumed that the 
furthest-thrown fragments were found and recorded. This is generally a 
reasonable assumption in reports on accidental explosions since the 
furthest-thrown fragments are generally of sufficient size to be observed 
and found, and since, in the thickly strewn close-in region, only large 
fragments or fragments considered significant in determining the cause of 


the explosion are recorded. 


On the general shapes of the actual and theoretical curves are 
similar; this similarity is quite pronounced for the Pantex explosion. The 
average curve of fragment distribution for the Wolf Creek and Kankakee ex- 
plosions is observed to fall between the theoretical and actual curves for 
walis and roofs -- a reasonable expectation. The steeper rise of the Pantex 
curve below a range of 0. 8-0. 9 times maximum may result from the nature 


of the structure. 


It would be desirable to use these curves along with the maximum 
debris distance curves presented earlier to develop areal density patterns 
or probability of hit. To do so would require that both (1) the number of 
fragments in the outermost circumferential sector and (2) the nature of the 
curve in the region where it starts rising rapidly, be defined better. These 
factors may be related to the size and strength of the structure and their 
determination would thus require extensive experimental investigation. The 
currently available data, considered above, could be used for this purpose 
although it should be recognized that they are far from well defined experi- 


mental input. 
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2.3.2 Debris Dispersion Patterns of HE Explosions 


Several prior studies of debris dispersion about accidental and 
experimental explosions have been made. One was performed by Colonel 
Clark S. Robinson who studied the results of six ordnance explosions, then 
compiled and compared their results (Ref. 12). This study is reviewed 
here along with additional plots of the recorded data which were prepared 
in an effort to find greater consistency of the plotted patterns. The second 
study was made at the Pantex Ordnance Plant at Amarillo, Texas, where 
about 31,000 fragments totaling 85,000 lb from a reinforced concrete struc~ 
ture were located and weighed after a planned interior explosion of 2,000 lb 
of high explosives (Ref. 11). The raw data available from this study were 
used in making a detailed study of the debris dispersion pattern from this 


explosion. 


Army-Navy Explosives Safety Board Investigation 


of Debris Dispersion 


Colonel Robinson plotted the specific area per missile against 


ground range. Specifically, the following explosions were studied: 


Badger Ordnance Works, 1945: 


7,500 lb of nitroglycerine exploded accidentally 
in a barricaded storehouse building of light frame 
construction. Data on about 600 fragments from 
building and contents were recorded. 


Cornhusker, Nebraska, 1945: 


10, 000 lb of explosives detonated accidentally in 
a loading plant manufacturing bombs. Chief 
interest in the investigation centered about 
fragments of machinery and contents and few 
building fragments were included. 


Portage, Ohio, 1943: 


Bombs amounting to 40,000 lb of HE detonated 
accidentally in an arch-type, earth-covered, 
igloo magazine. Data on whole bombs and 
bomb fragments constituted the bulk of recorded 
debris. 
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Umatilla, Oregon, 1944: 


64,000 lb of HE in bombs exploded accidentally 

in an Army-type igloo. Debris described con- 

sisted chiefly of pieces of concrete and frag- 

ments of machinery stored in the igloo. 
Hastings, Nebraska, 1944: 


100, 000 lb of HE ignited in a building where bombs 
were being loaded. Interest lay chiefly in large 
concrete fragments from the builidng itself, and 
only information on those weighing over 25 lb 
was included. 

Arco, Idaho, 1945: 


250, 000 1b of explosive in bombs was detonated 

in a planned explosion in an igloo-type magazine. 

Data on over 13,000 fragments were recorded. 

Specific area is plotted against ground range for these six 
explosions in Fig. 2.30. Several limitations of these plots, stemming 
from existent debris recording practices, must be recognized. First, the 
investigators made no attempt to record all fragments, especially in the 
first five explosions listed above; this tends to give higher values for 
specific area than actual, perhaps by several orders of magnitude. Secondly, 
in the thickly strewn close-in region, explosion investigators tend to record 
the exceptional fragments only; hence the reversal of several curves at 
short ground ranges in Fig. 2.30 is unrealistic. It is stated in the ASESB 
report, however, that at the greater ground ranges, it was customary for 
investigators to record all missiles of sizes that would damage structures 
or kill personnel -- that is, all fragments weighing one pound or more. It was 
observed that in these explosions the far-flying fragments were usually 


large and most were discovered and recorded. 


The individual curves of specific area in Fig. 2. 30 do exhibit 
a consistent and characteristic shape. The relative positions and crossing 
over of the four intermediate explosions (10,000-lb, 40, 000-lb, 65, 000-1b 
and 100, 000-1b curves) make it difficult to derive a general expression 
for the functional reltaionship between specific area and ground range for 


contained explosions. This stems from the fact that the curves are based 
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on incomplete fragment counts, and that the actual specific area is a 
function of the nature of the structure itself, as well as the explosive 


weight and ground range. 


Several additional plots of the data from Fig. 2. 30 have been 
prepared in an effort to find general curves for these functions with a better 
fit. These additional plots, shown in Fig. 2. 31 through 2. 36, were derived 
by selecting five points on each of the explosion curves of Fig. 2.30 and 


computing new expressions as follows: 


Fig. 2.31 First-Order Logarithmic Curves of Debris 
Dispersion for Selected Explosions. 


Log-log curves were fitted through the five selected 
points by the least-squares method, producing a 
series of curves of the form: 


logig An =k [ogy p}° 


No consistent relationship between either the slopes 
{n) or the intercepts (k) and the equivalent yield 
are apparent from these plots. 


Fig. 2.32 Second-Order Logarithmic Curves of Debris 
Dispersion for Selected Explosions. 


Log-log second-order curves were fitted through the 
same five points as above for each curve. Distinct 
separations between the curves for the various yields 
are absent, and the negative slopes near ground zero 
are Questionable. 


Fig. 2.33 Modified Second-Order Logarithmic Curves 
of Debris Dispersion for Selected Explosions. 


Log-log second-order curves were fitted, again by 
the least-squares method, through six points for each 
curve -- the five used previously and the arbitrary 
addition of point (1,1) as a data point. The resulting 
curves exhibit some separation at the shorter ground 
ranges, but the crossing-over of curves and their 
relative positions cannot be explained. 
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2.34 Parabolic Curves of Debris Dispersion 
for Selected Explosions. 


These curves are parabolas of the form 


2 

An = kD 
where, for each curve, k is the average value 
for a series of parabolas through point (0, 0) and 
each of the five selected data points. The slope 
appears adequate for some curves (notably the 
7, 500-1lb and 250, 000-1lb explosions), but the rela- 
tive positions of the curves cannot be explained 
in terms of equivalent yield alone. Total quantity 
of available material for fragmentation probably 
is involved. 


2.35 Logarithmic Parabolas of Debris Dispersion 
for Selected Explosions. 


These curves are of the form 
2 
Logig Ap =k (log) g D) 


where, for each curve, k is the average value for 
a series of parabolas through point (1,1) and each 
of the five points on the curves. 


2. 36 Second-Order Semi-Logarithmic Curves 
of Debris Dispersion for Selected Explosions. 
These curves are of the form 


Log) 9 An = atbD+t eb? 


where the coefficients a, b, andc are determined 
by the least-squares method. 


Of the six sets of curves, the logarithmic parabolas appear to 


provide the best fit. Relative positions of the curves cannot be explained, 


i.e.» why they should not be consecutive in the order of yield. As stated 


earlier this may be a function of the amount of material available for 


fragmentation and the degree of fragmentation, which can be expected 


to vary with the quantity of explosives or the impulse. Because of the lack 


of consistency in these results, no further analysis was made of these results. 
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10, 000, 000 


Plotted curves are semi- 
logarithmic parabolas of the form 


Log sat bD+ cD* 


100, 000 — lo D 
Where the coefficients a,b, and c 


are determined by the least- 


squared method 


Specific Area, (sq ft/fragment) 


ny 7,500 lb HE 
ao 10,000 lb HE 
e 40, 000 lb HE 
° 
° 


65, 000 lb HE 


100, 000 1b HE 
x 250, 000 lb HE 


0 1, 000 2, 000 3, 000 4, 000 5, 000 


Distance from Explosion, (ft) 


Figure 2.36 Second-Order Semi- Logarithmic Curves of Debris 


Dispersion for Selected Explosions 
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It will be recalled that investigators of accidental explosions 
tend to make incomplete counts of the smaller fragments in the thickly 
strewn close-in region about the source. Since plotted points in the fore+ 
going study were all taken at equal value, this would result in curves which 
are too shallow, and thus tend to overestimate the debris problem when 


extrapolations are made beyond the limits of the data. 


One measure of the fit of the fragment dispersion function can be 
obtained visually from data on the 1945 Badger Ordnance Works explosion of 
7,500 lb of nitroglycerine from the plot of data points as shown in Fig. 2. 37. 
Plotted points in Fig. 2.37 show debris dispersion expressed as specific 
area for each 10-ft interval in ground range from the point of burst. The 
average line drawn through these points was not originally computed by the 
least-squares method, but was an average based on the four central curves 
of Fig. 2.30 plotted at one-half the ground range. Though this original 
method of plotting was used merely because the Badger explosion was a 
barricaded structure, the line does appear to be a fair representation of 


this particular set of data. 


The two curves denoting the one-standard-error limits have been 
plotted visually to include two-thirds of all plotted data points between them. 
The log value of the standard error is thus measured as + 0.328, giving a 
standard error of ~-> 2.13 about the central line, assuming it were a true 


average and the distribution were normal. 
Debris Dispersion from a Reinforced Concrete Structure 


The 1960 planned explosion at Pantex Ordnance Plant provided 
extensive debris dispersion data (Ref. 11). This explosion involved the 
detonation of 2,000 lb of HE in the form of encased warheads placed inside 
the standard one-foot-thick reinforced concrete walls of a U-shaped bay. 
Two views of the structure are shown in Fig. 2. 38 and 2.39. The structure 
was completely destroyed. Debris was dispersed over a large area; the 
maximum debris distance being about 1,500 ft from the point of burst, 
Figure 2,40 is a post-shot view of the close-in region. The area was 
canvassed and all fragments found were listed by terminal location. About 


31,000 conerete fragments with a total weight of about 85,000 lb were 
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Figure 2.37 Debris Dispersion and Standard Error 
for 1945 Badger Ordnance Works Explosion 
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Figure 2.41 Area of Fragment Distribution 


from Pantex Ordnance Plant Test Structure 


95 


SECRET 


recorded along with their individual weights and terminal locations. Terminal 
locations were noted according to the serialized 50-ft square zones shown 

in the missile map, Fig. 2.41. As the missile map indicates, debris 
dispersion was non-uniform with the greatest concentration being to the 

south -- the direction of the front wall from the point of burst. Debris 


data were not analyzed in the original explosion report. 


Since the fragment list reported for the Pantex Ordnance Plant 
event constituted the most extensive compilation of data found on dispersion 
of fragments, these results were studied extensively to obtain data on 


fragment-size distributions and fragment dispersion. 


Over-all dispersion of fragments from the Pantex event is 
tabulated in Appendix D and plotted in Fig. 2.42 through 2.47. Two 
measures of dispersion are included here, specific area in sq ft per fragment, 
and areal dispersion in sq ft per lb of debris. The following plots have been 
prepared as alternate means of describing the dispersion function: 
Fig. 2.42 Log-Log Plot of Specific Area Vs Ground Range, 
with Second-Order Least-Squares Regression Line 


Fig. 2.43 Semi-Log Plot of Specific Area Vs Ground Range, 
with First-Order Least-Squares Regression Line 


Fig. 2.44 Semi-Log Plot of Specific Area Vs Ground Range, 
with Second-Order Least-Squares Regression Line 


Fig. 2.45 Log-Log Plot of Areal Dispersion Vs Ground Range, 
with Second-Order Least-Squares Regression Line 


Fig. 2.46 Semi-Log Plot of Areal Dispersion Vs Ground Range, 
with First-Order Least-Squares Regression Line 


Fig. 2.47 Semi-Log Plot of Areal Dispersion Vs Ground Range, 
with Second-Order Least-Squares Regression Line 
In each case the dispersion measures, (Specific Area or Areal 
Dispersion), have been computed for 50-ft and 100-ft wide concentric 
circular bands about the point of detonation. Computed points were plotted 
at the midpoints of the circular bands. Since the original debris data from 
this event were collected for 50-ft squares from coordinates near the point 
of burst, the approximation of the actual dispersion function is based on 
considering the debris within any 50-ft square to be within the circular band 


containing its center. Regression lines are computed on the basis of all 
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Figure 2.43 Debris Dispersion for Reinforced Concrete Structure 


Semilog Plot with First. Order Regression Line 
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Figure 2.45 Areal Dispersion for Reinforced Concrete Structure 


Logarithmic Plot with Second-Order Regression Line 
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points for the 50-ft circumferential bands as listed in Table D-2. 


Variation in average fragment weight with ground range is 
plotted in Fig. 2.48. Again, these values were computed for concentric 
circular bands of 50-ft and 100-ft widths. Though the largest fragment 
recorded weighed about 4,000 lb, it is seen from this chart that the average 
weight of fragments at virtually all ground ranges was less than 12 lb. 

The bimodal nature of this curve is unexplained, and since this computation 
includes concrete fragments only, it does not stem from a two-phase 

nature of material. Relative distances from explosive sources to the various 
structural panels and the shielding effects of adjacent bays to the east of 


the bay containing the charge may have influenced this. 


Fragment-size distribution for the 50-ft and 100-ft debris 
zones is plotted in Fig. 2.49 and 2.50. Figure 2.49 shows the total number 
of fragments above any indicated weight which were found at various ground 
ranges. This figure shows that relatively few fragments above three pounds 
were found at any ground range. An improved indication of variation in 
fragment-size distribution is obtained by plotting the cumulative percentage 
of fragments above indicated weights, as shown in Fig. 2.50. There is 
little separation between the cumulative-percentage curves over much of the 
range of fragment sizes in this figure. The pronounced ''flattening-out" of 
these curves at the 3-lb fragment weight suggests that an optimum design 
point for debris hazards may exist, providing of course that these curves 
are truly characteristic. At all ground ranges, less than 6 percent of total 
fragments was above three pounds in size. It is interesting to note that at 
the greatest ground ranges (1200-1500 ft) no fragments above one pound 
were found; but at the intermediate ground ranges (600-1200 ft) the percentage 
of fragments over three pounds in size was greater than at the shorter ground 
ranges (0-600 ft). The largest sized fragments (over 100 1b) were found only 
at the shorter ground ranges. The tendency of making incomplete counts of 
the smallest fragments in the close-in region would actually result in the 
cumulative percentagecurves of Fig. 2.50 being overstated for the greater 


fragment sizes. 


SECRET 


102 


| 
wa fH ete GE EEN 
| 


Espueq [eljUsTasUINIAIS 13 OG UO peseg _ 


fore ee AE a Lan nee eos f- mt 7 
5 yeiquarToyMIN IAD iW AY uo posed od [’ 
+ V7 x 1 


es 


| W Tae 


MEINEM TT 


NG 
| 
[ 
| 
| 
+ 
| 
| 
| 


(staqap jo qt/3 bs) ‘uotsradsiq teary 


000 ‘Oot 


000 ‘000 ‘I 


000 ‘000 ‘2 


Figure 2.46 Areal Dispersion for Reinforced Concrete Structure 


Semilog Plot with First-Order Regression Line 
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Figure 2.47 Areal Dispersion for Reinforced Concrete Structure 


105 


Semilog Plot with Second-Order Regression Line 
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Fragment dispersion for various weight classes is plotted in 
Fig. 2.51 and 2.52. Actual counts of fragments are presented in Fig. 2.51. 
it is seen here that the heaviest fragments did not travel to the extreme 
ranges. Fragments of 1,000 lb and heavier were found only within 200 ft, 
those of 100 1b or more were found only within 600 ft, and those of 
10 1b or more were not found beyond 1200 ft. This is contrary to what 
would normally be expected on consideration of air drag as shown in 
Chapter Five. The comparison is not altogether appropriate, however, 
since in this actual structure, the fragments which are subject to 
forces of sufficient magnitude to cause large motion are also likely to be 
broken up to a greater extent. Proportionate distribution of fragments of 
various size classes, among the various ground ranges, is shown in 
Fig. 2.52. Other than for the extremely heavy fragments (over 100 lb) 
the proportionate distribution did not change much for the different 


fragment weights. 
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Figure 2.51 Fragment Distances for Various Weight Classes 


for a Reinforced Concrete Structure 
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Figure 2.52 Fragment.Distance Distribution for Various Weight Classes 


for a Reinforced Concrete Structure 
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CHAPTER THREE 
FRAGMENTATION, EXPERIMENTAL OBSERVATIONS 


To predict the hazards of debris from nuclear explosions, in 
addition to defining the maximum range of debris and the distribution of 
debris within this limit, it is well to define the nature (i.e., size or weight) 
of the expected missiles, Fragment weight combined with fragment velocity 
determines energy level of the fragment, and thus the loading upon equip- 
ment of personnel struck by the fragment. No definitive experimental 
investigation of structural fragmentation, relating fragment-size distribution 
to structural strength and loading parameters has yet been made. Recourse 
was made to collecting and summarizing past experimental investigations 
in fragmentation and to describing such general behavior and characteristic 
patterns as have been developed. Five such investigations are described 
in this chapter; research of the British Coal Utilization Research Association 
on coal breakage from random forces, research of the Safety in Mines 
Research Establishment of Great Britain on explosively detonated stone 
blocks, Stanford Research Institute model tests of containment structures 
fractured by internal explosions, an extensive study of fragments froma 
planned explosion of an ordnance structure at Pantex Ordnance Plant and 
findings of Project 4. 5 of Operation JANGLE. 


The fragmentation of materials produced by mine charges or 
ore crushers has been the subject of considerable study. Many attempts 
have been made to analyze the fractions produced by these processes. It 
has been shown experimentally that the higher the loading on the source 
material, the smaller the fragments produced, and that a wide range of 
fragment sizes are produced at any loading. An "Ideal Law of Breakage" 
which shows excellent fit to experimental data has been developed from 


coal-crushing investigations. 


The extensive collection of data on the sizes of more than 
30, 000 concrete fragments from the Pantex Ordnance Plant event afforded 
the opportunity for a detailed analysis of fragment-size distribution. An 


interesting result of this particular study is that only about 3 percent of 
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all the fragments produced weighed more than three pounds and that these 
fragments accounted for nearly 75 percent of total weight of all fragments 
recovered, Thus, as structures or equipment are built with increased 
resistance to fragments, the probability of being hit by a fragment 


sufficient to cause damage declines substantially. 


The only structural fragmentation study conducted on full-scale 
nuclear tests was limited work performed on reinforced-concrete wall 
panels erected over the crater zone in Project 4.5 of Operation JANGLE 
(Ref. 5). Size distributions of the larger fractions were plotted as part of 
this project and it was noted that the JANGLE data did not preclude the 
possibility that the fragment-size distribution of concrete source material 
caused by the underground nuclear shot followed the same pattern as coal in 


a mine or ore ina crusher. 


3.1 Fragmentation of Coal 


A number of experimental investigations of coal fragmentation 
characteristics have been made. The mining industry, interested in mini- 
mizing dust formation and in producing fragments of an appropriate size 
for handling, has conducted numerous investigations of fragmentation from 
blasting and crushing operations. Assuming a brittle material with a random 
distribution of internal weaknesses, it has been deduced that when a single 
lump is fragmented by forces sufficiently violent to make breakage equally 
likely at any point, the broken product should have a fragment-size 


distribution obeying the exponential law (Ref. 13): 


M(x) = 1 -e7*/* 


where 


M(x) = percentage of material smaller in size than x 


x/x dimensionless measure of fragment size, such as 


ratio of a characteristic length to a mean length. 
This has been called the Ideal Law of Breakage. Experiments with large 
lumps of coal broken under conditions approximating random fracture tend 
to conform to the exponential law, at least for the smaller fragment sizes 


(under about )/2-in. equivalent spherical diameter). Fitting data to the 
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equations 


M(x) = c¢(1 - eX! *) 


and = 
M(x) = 1- el x, 


researchers of the British Coal Utilization Board plotted experimental frag- 


mentation data on broken coal of various types as shown in Fig. 3.1 (Ref. 13). 


3.2 Fragmentation of Explosively Detonated Stone Blocks 


In research conducted by personnel of the Safety in Mines 
Research Establishment, stemmed charges of l~oz and 2-0z of coal- 
mining explosive were fired in stone blocks 18-in. in diameter by 30-1n. 
long contained ina steel chamber (Ref. 14). Weights of the full-size- 


distribution fractions, including all dust, were measured. 


Weights of the various fractions have been recomputed in terms 
of cumulative fragment quantities in Appendix E and are plotted in Fig. 3. 2. 
The cumulative curves all follow the same general trend, which 1s nearly 
linear on logarithmic coordinates. Slopes of the lines vary with the quantity 
of charge, and seemingly with the diameter of the shothole, both of which 
would be determinants of the impulse impinging on the stone blocks. Average 
slopes of the cumulative size-distribution curves for 2-oz charges were all 
significantly greater than those for l-cz charges and showed from only 
one-fifth to one-tenth the proportion of larger fractions. Thus, the larger 
quantities and smaller sizes of fragments were produced by detonations of 
the larger charges providing greater impulses -- certainly not an unexpected 
result. This tends to lend support to theories, that, in impulsive frag- 
mentation, new surface area produced through fragmentation may be 


quantitatively related to explosive impulse or energy. 


3.3 Fragmentation of Concrete Shielding of Reactor Models 


The Stanford Research Institute conducted experimentation on 


model concrete biological shielding poured directly in contact with model 
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pressure vessels (Ref. 15). The models, shown in Fig. 3.3, consisted 
of a reactor vessel, the concrete biological shield, and a stepped plug 
filling the access opening at the top of the vessel. Vessels were filled 
with varying amounts of water (75 and 100 percent) and subjected to the 
detonation of scaled energy sources centered radially and axially in the 
pressure vessel. The following model tests were selected for detailed 
plotting in this study, since individual fragment weights were available 
for these events. Original source data from these tests are included 


in Appendix F. 


Table 3.1 


MODEL TESTS SELECTED FOR PLOTTING 


: Water Number of 
Model Energy in Fragments 
Test Megawatt- | Source Vessels, | Above | Tota 
No. seconds Material (%) 0. 1 lb | Recorded 
15 Pyracore 1 75 14 24 
17 Pyracore 1 100 7 9 
14 Pyracore 1 75 2 2 
i 16 MDF 1 100 41 373 
19 MDF 1 75 | 55 287 
} 
5 Pyracore 1 160 { 23 63 
3 Pyracore 1 75 45 96 


The cumulative total number of fragments above the indicated 
size is plotted in Fig. 3.4 and 3.5. Figure 3.4 includes all recorded 
fragments, while Fig. 3.5 includes only those fragments weighing more 
than 0.1 lb. Three explosion factors in the test setup which influence the 
loading are the percent of water in the vessel, the quantity of explosive 
(as indicated by the magnitudes of the simulated event), and the equivalence 
of the explosive source material. Tests, numbered 16 and 19, in which 


MDF was used as the explosive source material produced the most 
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Concrete Shield 
a Maximum Diameter 12 in. 
Height 20 in. 


Steel Plug 

Maximum Diameter 5-1/2 in. 
_ Minimum Diameter 4-1 /2in. 
/“ Height 4in. 
Weight 221b 


Section A-A 


a Screw Top 
_ 4in. OD 
1/4in. Thick 
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1/2 in. x 1/4 in. 
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Figure 3.3 Model Reactor with Concrete Shielding 
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Figure 3.4 Cumulative Fragment Distributions for 1/24-Sca.e Shielded 
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distribution tabulated in Appendix D. 
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complete fragmentation. This certainly appears to be a function of the 
impulsive loading (or energy input) since the amount of gaseous reaction 
produce in the energy source is five times as high for MDF as for 
Pyracore (100% vs. 20%). 


Other than for Test 15, it was observed that fragmentation 
became more complete with increasing magnitude of the simulated event, 
higher water level in the vessel, and increased quantity of gaseous 
combustion products from the charge. These are all factors which would 
increase the loading on the concrete shielding. ‘In Test No. 15, 24 frag- 
ments were produced and 85 percent of the total concrete weight was in 


the largest fragment. 


A highly favorable result of this series of tests is the consistency 
in the shape of the cumulative curves of Fig. 3.4 and 3.5. Crossing-over 
of curves exists 1n these plots, but the general shape of the curves is 
consistent and the relative position of the successive curves appears to 


bear some relation to the impulsive loading. 


Similar results were observed ina series of 1/12-scale shielded 
reactor model tests (Ref. 16). Only the magnitude of the event simulated 
was varied in this series of tests, (data are tabulated in Appendix G and 
plotted in Fig. 3.6). As Fig. 3.6 shows, separation between the cumulative 
distribution curves for different events is decidedly pronounced, and the 
general shape of the curves is quite consistent. This series of tests appears 
to support contentions that fragmentation becomes more complete with 
increased loading and implies that the relationship between fragmentation 


and loading may be quantified, at least for simple ideal structures, 


Stanford Research Institute selected three criteria to assess 
the fragmentation and fragment dispersion in the reactor shielding model 
tests. The quantity of new surface area created in the fragmentation 
process was regarded as an indicator of the energy absorbed in breaking 
the concrete shielding. It can be shown that this newly created surface 
area (termed fragmentation in the SRI reports) can be approximated from 


the individual fragment weights: 
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Fragmentation @ New Surface Area 


Pi ee 


where We is the individual fragment weight. 


Similarly, the throw, defined as the summation of the products of weight 
times distance thrown for individual fragments, was taken as an indicator 


of the total amount of input energy diverted into transport of fragments; 


Throw @ 2 WD, 


where D; is distance thrown for the individual fragment. 


A measure of structural integrity representative of the degree to which 
concrete shielding remained intact was defined as the ratio of the sum of 
the squares of individual fragment weights to the square of the total shielding 


weight; 


Le 
(= W) 


The inverse of this integrity ratio was considered to be a measure of the 


total damage to the shielding structure. 


The measures of fragmentation, throw, and integrity were 
computed for the various reactor shielding models by SRI. Plots of these 
values are included in Fig. 3.7 through 3.10. For the various series 
of models, the selected debris parameters (fragmentation, throw, and the 
inverse of the integrity ratio) show varying degrees of consistency in the 
shape of the plotted relationships. Some of the more pronounced deviations 
between plotted points and the curves as drawn are perhaps explanable in 


terms of the considerable scatter customarily experienced in explosion testing. 


Figure 3.7 shows results of 1/12- and 1/24-scale model tests 
for the debris criteria. In this figure the fragmentation and throw for the 
1/12-scale models are reduced by the reciprocal of appropriate scale 


factors: 
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Fragmentation: Awa » by 4 
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eee WINS ace sy 
Throw: W,Di~rA” Wi Di, by 8Y2 


where A is the ratio of scale factors. 


Comparison of debris parameters for model reactor shielding 
with different quantities of contained water are shown in Fig. 3.8. Some, 
but not complete, separation of data points is apparent in the series of 
model tests in Fig. 3.9 and 3.10. With further tests, more conclusive 


definition of these functions might evolve. 


No relationships can be established between the results of 
these reactor shielding model tests, fragmented by the impulsive loading 
of an internal explosion, and the fragmention behavior of other structures 
under nuclear blast loading. Under the longer duration nuclear loading, 
it is likely that peak overpressure may exert a greater significance than 
impulse on the fragmentation process. It is significantly demonstrated, 
however, that fragmentation patterns can be measured and related to loading, 
and that, for a uniform model structure consistent relationships are 
obtained. This suggests the feasibility of experimenting with suitable 
ideal structural elements, ona full-scale or model basis, to develop 


basic input data for debris hazard estimating purposes. 


3.4 Fragmentation of a Reinforced Concrete Ordnance Structure 


Fragment counts from the Pantex Ordnance Plant event described 
previously provided data on the individual weights of more than 31,000 
recovered concrete fragments with a total weight of more than 85,000 lb 
(Ref. 11). Individual fragment weights recorded ranged from 1/16 lb 
to 4,000 lb. 


The fragment-size distribution of the recorded fragments is 
shown in Fig. 3.11. Fragment classes in this figure are in a geometric 
progression, each class being twice as large as the preceding class. 
From this chart it appears that the greatest portion of fragments were 


in the weight range of 0. 250 to 0. 499 lb, with equivalent diameters of 
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1.8 to 2.2 in. The expected preponderance of relatively small fragments, 
under 4 lb or 4,4-in. diameter, is apparent. The relatively few fragments 
in the 0. 060-0. 249 lb range is questionable and could result from the 
tendency to not record all the smallest items at close-in ranges (see 

Fig. 2.51). 


Though the smaller fragments are produced in the greatest 
quantities, the amount of material involved is relatively small. The 
distribution of total weight of fragments, pictured in Fig. 3.12, shows 
more than 65 percent of the total material fragmented into pieces weighing 


between 32 and 255 lb, that is, with equivalent diameters from 8.9 to 17. 8 in. 


Cumulative distributions of fragment sizes are approximated in 
Fig. 3.13. The median fragment weight is found to be about 0. 22 lb, and 
more significantly -- more than 95 percent of all fragments weigh about 
3 lb or less, i.e., have equivalent diameters of 4 in. or less. Less than 
1 percent of all fragments weight 90 lb or had equivalent diameters greater 
than about 1.1 ft. In general, fragments above one pound in size are 
considered potentially lethal to personnel. Only about 20 percent of all 


recorded fragments from this test were above this size. 


Cumulative weight distribution for fragments is presented in 
Fig. 3.14. Whereas the greatest number of fragments produced in this 
event were about 3 lb or less in size, Fig. 3.14 shows that somewhat over 
50 percent of the total weight of fragments was in boulders weighing greater 
than 70 lb or having equivalent diameter of about 11.5 in. or more. Over 
70 percent of total fragment weight was in pieces above 3 lb in size. Thus, 
although substantial quantities of concrete fragmented into large pieces, the 
number of these was small. The probability of being bit by small pieces 


is much higher than that of being hit by large boulders. 


3.5 Fragmentation of Concrete Walls 


The only fragmentation study of structural elements under 
nuclear loading was conducted as part of Project 4.5 of Operation 


JANGLE. For this event a series of six reinforced concrete wall panels was 
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erected over the crater zone of the underground shot as shown in Fig. 3.15. 
This event involved a 1. 2-KT device emplaced at a depth of burst of 17 ft; 
the apparent crater diameter was about 260 ft. The wall panels were thus 
placed fairly close-in relative to the total crater -- i.e., at ground ranges 


from 14 percent to 40 percent of apparent crater radius. 


Blue 


Yellow 


> an eae 


Charge Centex 


Figure 3.15 Placement of Wall Targets in JANGLE U Event_ 
The fragment-size distribution of wall fragments reported 

from this test is presented in Fig. 3.16. These curves are plotted 
through the data points at a slope of 0. 5 (the slope predicted on the basis 
of the Ideal Breakage Law). The JANGLE report pointed out that only the 
larger fragments were collected and that these would normally be expected 
to deviate from the ideal fragment-size distribution when a limited 
amount of material is available. It was concluded that the JANGLE data 
did not preclude the possibility that the breakup of concrete source material 
caused by a nuclear detonation follows the same pattern as coal in a mine 


or ore ina crusher. 
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It is noted that the relative positions on the curves (and their 
associated data points) are not arrayed in the same sequence as the 
initial ground ranges of the walls. The lowest curve, representing the 
largest fragments, corresponds to the wall closest to surface zero, where 
vertical components of the load may be expected to be higher. However 
the two topmost curves, representing the smallest fragment sizes, correspond 


to the middle two walls rather than the farthest-out. 
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CHAPTER FOUR 


FRAGMENTATION, ANALYTICAL CONSIDERATIONS 


4.1 State of Knowledge of Fracture Mechanics 


Analytical consideration of debris resulting from the fracture 
of concrete or masonry structures poses two fundamental questions: 1) What 
is the fracture strength, i.e., the ultimate ] ad which the structure can 
sustain? and 2) What is the number and size of the resulting fragments ? 
The first question is by far the easier to answer and it is discussed briefly 


from the point of view of existing fracture theory. 


The theoretical strength of materials is of the order of 100 to 
1000 times the observed strength. Griffith (Ref. 17) proposed that this 
difference can be rationalized in terms of pre-existing flaws contained 
in the solid. This model of a solid containing an array of flaws is the basis 
of the "fracture mechanics" approach to the problem. Treatment of this 
problem is concerned with the growth of pre-existing flaws contained in the 
solid. This model of a solid containing an array of flaws is the basis of the 
"fracture mechanics'' approach to the problem. Treatment of this problem 
is concerned with the growth of pre-existing flaws and conditions of 
instability which can change the crack propagation from a slow process to 


that characteristic of a fast-running crack. 


Griffith formulated the condition of stability, under load, ofa 
body containing a certain type flaw upon two theorems. The theorem of 
minimum energy states that the equilibrium state of an elastic body, 
deformed by specific surface forces, is such that the potential energy of 
the whole system is a minimum. He obtained his new criterion of rupture 
by adding to the theorem of minimum energy the statement that the equil- 
ibrium position, if equilibrium is possible, must be one in which rupture 
of the solid has occurred, if the system can pass from the unbroken to 
the broken condition by a process involving a continuous decrease in 
potential energy. Thus, in the propagation of a crack, stored potential 
energy is released, but the potential energy of the system is increased 


by the creation of new surfaces (surface energy). Griffith's condition for 
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continuing propagation of a crack is that the resultant potential energy of the 
system is decreasing. Likewise, the equilibrium crack size is one in which 
the decrease in potential energy just equals the increase in surface energy. 
The theory predicts strength reasonably well for bodies such as glass and 
ceramic which behave in a brittle fashion. Using the Griffith theory, we 
can reason that there will be a distribution of strengths in a given specimen 
in the sense that a different amount of force will be needed to fracture a 
specimen at one point than at another. If one assumes that the flaws are 
distributed at random with a certain density per unit volume, then 

the statistical formulation of the strength problem becomes apparent. The 
strength of a specimen is determined by the weakest point in the specimen 
or by the smallest value to be found in a sample of size n where n is the 
number of flaws. Clearly, n increases as the volume increases and, there- 
fore, the problem of finding out how the strength depends on the volume of 
the specimen is equivalent statistically to studying the distribution of the 
smallest value as a function of n, the sample size. This statistical prob- 
lem is an important one on which much theoretical work has been done. For 


breaking strength the major contributions have been made by Weibull (Ref. 18). 


Essentially, Weibull assumes that the probability of failure of a 


unit volume of material can be represented by a distribution function of the 


form 
m 
| o-0 
F(c),= 1- ex - dv o 20 4.1 
(6), P AZ 20, (4.1) 
= 0 o <o 

where u 

F(c) = failure probability for stress o 

Tem = constants of the material 
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Once the unit probability of failure F(c) is known, it is a reasonably 
straightforward procedure to find the probability of failure of any 
structure under any known system of stresses. Note, however, that the 
failure mode is not unique since the location of the weakest link is a 
statistical quantity. Since the likelihood of failure is greatest in regions 
of high stress, it is certainly possible to anticipate the origin of failure. 
This predicability is relied upon when brittle test specimens are designed 
to break in the gage length. But even here it is quite common to get 


fractures outside of this region. 


The applicability of the weakest link theories, and especially ‘. 
the Weibull theory, to the prediction of fracture strength in masonry 
and ceramics is a problem which is currently undergoind the most intensive | 
investigation. Since the theories, at best, only treat the static load problem, 
their use in the exceedingly complex debris formation problem is simply 


too much to expect. 
4.2 Mathematical Model 


Since a sophisticated treatment of fragmentation involves 
application of fracture mechanics to an extent that is well beyond the current 
state-of-the-art, our objective here is to formulate a relatively simple 
mathematical model which can be used to predict debris formation and 


which can be adapted to include advances in fracture mechanics. 


The response of structures to nuclear blast loading is calculated 
by considering the response of an equivalent single-degree-of-freedom 
system (equivalent to neglecting all but the fundamental mode of vibration 
of the structure). The stiffness of this equivalent system is assumed to be 
elastic - perfectly plastic so that the structure's response into the plastic 


range can be considered. 


Since this model has proved adequate for analyzing the response 
of structures up to failure, we propose a simple extension of the model 
to include fracture. The model is shown on Fig. 4.1. A structure is 


represented as an equivalent mass M, and resistance R(x). The 
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F(t) 


M. 
x 
R(x) 
(a) Model 
R(x) F(t) 

F 

o 
f Fracture 

ae X¢ x a t 
(b) Resistance (c) Forcing Function 


Figure 4.1] Fragmentation Model 
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resistance is taken to be elastic-perfectly plastic up to a displacement (x,) 
at which the structure fractures. Such a model can be used to predict 
the velocity and time at which the equivalent mass reaches the fracture 


displacement. 


Limitations to this model are obvious. First, the mass distri- 
bution of fragments from the structure will not be determined by this model. 
Secondly, it is somewhat presumptuous to assume that all fragments are 
formed at a single displacement. These are questions however, which must 
be answered by advances (experimental and analytical) in fracture mechanics. 


Here, we say 


M, = aM 


f e 
where, 
M, = mass distribution of fragments 
a = normalized mass distribution of fragments to be 
determined. 


Similarly, there would be a statistical distribution about Xr 
which would result in a corresponding distribution in the time of fracture 
and velocity spectra of the fragments. Thus, refinements of the basic 
simple model would be statistical in nature. An experimental program 
could be undertaken to establish the nature of these distributions for 
particular types of structures (e.g., brick walls, concrete structures, steel 
frames). Hence, it is feasible to obtain answers to these questions (at least 
insofar as required to modify our model) without waiting for a complete 


understanding of the fracture mechanics. 


The equation of motion for the model of Fig. 4. 1(a) is, 


2 
M. es + Rx = F(t) (4. 2) 
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For convenience, the following nondimensional parameters are 


defined: 
T 2 t/t, 
oc = x/X, 
a2 2c oe ft)? 
e*o - 
w = circular frequency of elastic system 
2 F 
GO ee 
R 
x2 
= f 
Pp" 
e 


The equation of motion then becomes, 


on +p? g= p87 (1 - 7) 


for g<1 (4. 3) 
and 

on + Be = B* 87 (1 - 7) (4. 4) 
for 

1<G Sp 


where the double prime denotes differentiation with respect to T. 


Equation 4. 4 has a solution of the form (for zero initial conditions), 
t= 8? [-cos pr + RET ties] (4. 5) 


for 
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and 
Sane 2,2 
_ §° (6° - 1) 2 pv 2 2 3 
ig eee CT et) ? tT - 34 (T,) ae (4. 6) 
1 
+ o (Tt - T) +1 
for 
1<E< p 

where 

1 = time at which © = 1; determined from Eq. (4. 5). 

Bt = nondimensional velocity at TQ) determined from Eq. (4. 5). 


Solutions for Eq. (4. 5) and (4. 6) were obtained for final velocity 
and the time of fracture. These solutions are presented graphically in 
Fig. 4.2 through 4.4. Thus, the velocity of fragments and time of fracture 
can be found from this analysis. These would then be inputs to the analysis 
of debris transport wkich would provide a complete displacement history 
of the fragment. 


4.3 Example 


As an example of the application of this analysis consider the 
wooden siding from a railroad car. The failure of such siding was observed 
during the UPSHOT-KNOTHOLE series of weapon effects tests (Ref. 19). The 
velocity at failure of siding located in the 10-psi overpressure region was 


observed to be 60 fps. 
The ultimate elastic resistance of a simply supported rectangular 


beam is, 


Bins % 2 b (1. 25) x 1400 


o 7 L y x12 rae 
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Figure 4.2 Motion of Mass At Failure (6 = 5) 


SECRET 


145 


SECRET 


Ax1d0T 3A, 


of Mass At Failure (§ =10) 


Figure 4.3 Motion 
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Velocity 
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Frequency, (t ) 
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Figure 4.4 Motion of Mass At Failure (6 = 100) 
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where 


width of member, L 


span, and 


a 
i 


depth of member R total uniform load on beam. 


¢ = yield strength of material, 


The deflection of such a member is, 


3 
Roe _ 37x 10° x 1728 


x = = = 1 52 in 
REL © ee dee ir ee 
7 6.4 Ebd 64 x 10° x (1. 25) 


The circular frequency of a simply supported rectangular beam 


is 
2 { j 2 
_ oT E dg 
2 6 
= 100 | SS = 31.8 rad/sec 
where 


E = Young's Modulus 
= gravitational constant 


= density of material 


The 10-psi overpressure region corresponds (Ref. 20) toa 
2-psi dynamic pressure and a positive phase duration of 0.7 sec. The 


drag loading acting on the siding is then 


Fo =bL(2)cd=bx10x12x2x1= 240b 


Therefore 


and, 
B = wt, = 31.8x0.7= 22.3 
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Then if we interpolate between Fig. 4. 2 and 4. 3 and assume 


a fracture displacement of p= 20, 


c' = 300 


0.125 


Tf 


or, in dimensional form, 


x 
i e 1,52 x 300 
aD AS eee a, ae 
fo) 
t = 0.125x0.7= 0.09 sec. 


Obviously, this fracture velocity is dependent on the value assigned to op. 
Therefore, the reasonably accurate prediction of this velocity does not 
establish the validity of our model. It does show, however, that reasonable 


values of p lead to an acceptable prediction. 


The subsequent motion of the debris fragments is studied in 
Chapter Five. The significance of debris velocity at fragmentation Xe and 


time of fragmentation t, are also considered in subsequent chapters. 


f 
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CHAPTER FIVE 
DEBRIS TRANSPORT BY BLAST WINDS 


The motion of a particle acted on by the nuclear blast wind is 
considered with the fracture velocity and time of fracture forming the 
initial conditions of the transport problem. The model used here assumes 
that the force acting on the particle is proportional to the square of the 
relative velocity between the particle and air. The k.ast parameters are 
assumed to be constant over the range of travel of the debris, and it is 
further assumed that the effective lengthening of positive phase duration 
(i.e. , the so-called time correction) which is due to the debris motion 
in the direction of shock propagation can be handled by a simple adjustment 
of positive phase duration. These assumptions reduce the problem to the 
solution of a one-parameter nondimensional differential equation. Thus, 
much can be learned about the general behavior of flying debris. Without 
these two assumptions, it would be necessary to treat each weapon yield 
and placement as a separate problem and no general observations could 


be made. 


The general equations are derived and numerical solutions for 
many problems of interest are obtained. The application of these results 
to specific debris problems is outlined and general observations regarding 


the behavior of flying debris are detailed. 
5.1 General Treatment of Problem 


The motion of a piece of debris of arbitrary shape is first 
considered and then the problem is reduced to the point where the vertical 


and horizontal motions are uncoupled. 


Consider the motion of the body shown in Fig. 5.1. By equil- 


ibrium, 
mx = f,, 
my = f. (5. 1) 
6 = fy - x) +f (yg - xy) 
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where 
m = mass of body 
j = product of inertia of body. 
é@ = coordinate describing orientation of body with 


respect to x-y coordinate system 


f Centroid (x, Y.) 


Center of Pressure (x), yy) 


< 


Figure 5.1 General Debris Particle 


The forces (£9 fy) are given by 


a Car(s) 


fale p(u- x2 (5. 2) 


x 


ac 
-}! d(y) 
ff. ek So 


*2 
~ +m 
y py g 
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where, 


a = projected area in x and y direction 


CaGaqy drag coefficient in x and y direction 


p = mass density of air 
u = particle wind velocity 
g = gravitational constant 


Equations (5.1) are coupled, in general, by the dependence of 
acy on¢. Consideration of this dependency on ¢ makes exact descriptions 
of debris particles necessary. Weare, however, concerned with the gross 
behavior of groups of debris particles of many different sizes and shapes 
rather than a detailed treatment of specific shapes. Thus we restrict 
our attention to those particle shapes where aCy does not depend on ¢ 
and, hence, the equations of motion Eq. (5.1) ure uncoupled. The trajectory 
dispersion which is due to ¢-variation of aC, can be investigated separately, 
thereby leading to a dispersion function which can be superposed on the re- 
sults determined for the simpler problem considered here. Thus we 
consider only the first of Eq. (5.1). The following nondimensional parameters 


are defined: 


ones at 
= t 
te) 
aC 
Os Ut. g 
o°o m 
aC 
b= pS x (5. 3) 
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where 
Us = characteristic air particle velocity such that 
U = UU; hit) 
C= characteristic time of air particle velocity-duration 
curve 


Substituting Eq. (5.3) into Eq. (5.1) we have 


ek [onin-e]? (5.4) 


where the prime denotes differentiation with respect to T. 
Equation (5. 4) is a Riccati-type nonlinear differential equation which can be 


linearized by the substitution, 


to s' 
bis -k= (5. 5) 


Equation (5. 5) reduces Eq. (5.4) to 


2 


2 h’s 
=37 =O (5. 6) 


s'+k@hs'+86 


However, closed-form solutions can be obtained to Eq. (5.6) only for special 
classes of the forcing function h(t) (e.g., h(t) = constant leads to second- 
order differential equations with constant coefficients; h(t) = h/t leads to 

a Cauchy-Euler equation). The blast-induced wind is unfortunately not 

one of these forms. Therefore, we are forced to integrate Eq. (5. 4) 
numerically. The case of h= 0 is of interest and the solution for this case 


can be readily derived from Eq. (5. 5) and (5. 6) to be, 


C= -k log (t+ C,)+C, 
and note that, 


ci = oar « eae (5. 7) 
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where Cc): C, are constants of integration to be determined from initial 
conditions, 

Observe from Eq, (5.7) that the velocity only asymptotically goes to zero. 
This can be explained by recognizing that when the surrounding air is 
motionless the retarding force acting on the particle is proportional to the 
square of the particle velocity. Thus, solutions based on our model will 
predict infinite displacement at an infinite time. It will be necessary to 
place a maximum time limitation on the debris particles flight based on the 


time required for it to hit the ground. 


Numerical solutions to Eq. (5.4) are obtained for a wide class 
of problems of interest. Specifically, results are obtained for a wide range 
of aerodynamic coefficients, various initial conditions, and a few possible 


negative phase representations of the blast induced winds. 


We must consider first the form of the wind loading h(t), however. 


The particle velocity can be related to the dynamic pressure q by; 


1 2 
GES ape Pee (5. 8) 


During the positive phase of the loading, the dynamic pressure is 
q = 4G eet a sanyo (5. 9) 
When Eq. (5.9) is substituted into (5. 8) 
2q, 


us 5 e’ (1-7). (5. 10) 


Thus, in terms of Eq. (5.4), 


uz uy bh (7) 
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where, for 0<T<1 


q 
0 p 


2 (5.11) 


h(t) = e7 (1 - 7). 


u 


yam 
£q 
Variation of geo and to with weapon yield and distance from ground 


zero is shown on Fig. 5. 2 and 5. 3. 

Relatively little data are available regarding the negutive phase 
dynamic pressure (i. e., wind blowing toward ground zero). We have taken 
the negative phase dynamic pressure to be of the form, 


(t - 1) 


q = 4, sint (5. 12) 


so that wind ceases at t= 1 +k. 


It can be argued that, to have ambient conditions at ground zero 
some time after detonation of the weapon, the area under the positive phase 
dynamic pressure curve must equal the area under the negative phase curve. 


This results in, 


Io 
K = 0.34 ——. (5. 13) 
I 


The wave form descriptions apply to a fixed location with 
reference to ground zero. The debris particles, however, move with 
respect to ground zero. Note that uy and o. vary with distance from 
ground zero so that, to exactly represent the forcing function, the particle 
and pressure wave motion must be followed. As an example, consider 
a particle originally located at (x) and assume the pressure wave arrives 
at this location at time t= 0. The peak particle velocity and positive phase 
duration are uot) and t (x). At some later time (At) the particle has 
moved to (x + Ax); the pressure wave arrives at this location at some time 
less than (At), say (At'). Thus, the particle velocity at this time and at the 


actual debris location is given by, 
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u = u_ exp |- (At! - At)/t, .-| 
° 


“At! - At 
aay as (5. 14) 


where 


uy ts are evaluated at (x + Ax). 


Equations (5. 14) cannot be used to represent the particle velocity unless 
attention is restricted to specific weapon yields. Thus, to retain generality 
in the study, this refinement is not considered here. Note, however, that 
the form of Eq. (5.14) is well suited for inclusion in a numerical integration 


routine such as will be used to integrate the equations of motion. 
5.2 Numerical Results 


A computer program was written for the UNIVAC 1105 
digital computer to numerically integrate Eq. (5.4). The Runge-Kutte-Gill 
numerical integration procedure was used. ‘Solutions were obtained for 
zero initial conditions of the debris particle, and a complete range of 
realistic aerodynamic coefficients. The effect of negative phase wind on 
particle motion is also considered. It was the objective of this numerical 
analysis to obtain results in a form such that mass-velocity curves can be 
constructed for a given attach condition and location from a source of 


debris. 
5.2.1 Zero Initial Conditions Excluding a Negative Phase 


We consider first the piece of debris that is free to move, acted 
upon only by the positive phase wind loading. Solutions are obtained for 
values of © inthe range 0.1<9 <10. These are presented in Fig. 5.4 
and 5.5 in the form of curves of nondimensional velocity versus noridimen- 
sional distance from original debris source. Time is included as a 


parameter on these curves. 


As stated, attention has been limited to the horizontal motion 
of a particle. These solutions are invalid if the particle hits the ground. 
Therefore, the time parameter on the curves of Fig. 5.4 and 5.5 can be 


used to establish the range of validity of the results based on initial height 
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h of the debris and any vertical component of velocity Yo imparted to the 


particle. Based on elementary kinematics the particle will hit the ground 


at 
Yo 2h 
t = —! 1+ +] ; 5.15 
7 =“ ( ) 
fe) 


An interesting observation regarding the size and shape of 
particle that will travel the furthest can be drawn from these data. To this 
end, the parameter @ is plotted as a function of distance for various times 
on Fig. 5.6. Note that the results of this analysis predict that the size and 
shape (i.e., aerodynamic coefficient) of the piece of debris that travels 
furthest depends on the length of time required for the debris to hit the 
ground surface. The longer this time, the larger the piece of debris that 
travels furthest. This could explain the apparent randomness in high- 


explosive debris data. 


5.2.2 Zero Initial Conditions Including a Negative Phase 


As has been mentioned, little data which quantitatively describes 
the negative phase wind are available. A sine wave is assumed for the 
wave shape of the negative phase wind, and the area under the sine wave is 
made to equal the area under the positive phase wind. In this study, the 
peak negative phase wind velocity was taken to be 0.02 times the peak 


positive phase wind velocity. 


Velocity-distance curves for this negative phase wind are given 
in Fig. 5.7 and 5.8. The sharp departure of these results from those for 
zero negative phase is of considerable interest. Given sufficient time, 
all particles move toward ground zero. Of course, these required times 
of flight are sufficiently long that this reversal of velocity can probably 
not be realized for most practical cases. There have been instances, 
however, at the Nevada test site where pieces of structural debris were 
found closer to ground zero than the place they started. Of course, 
structural debris requires some time to be ripped from the structure so that 
a portion of the positive-phase wind impulse is not effective on the piece 
of debris. Thus the velocity reversal occurs sooner and can be realized 


with practical times of flight. 
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5.2.3 Effect of Initial Conditions 


The results described in the previous two sections must be 
modified to include the initial condition of failure time and velocity at failure. 
One could accomplish this modification by constructing velocity-distance 
curves (similar to Fig. 5.4 and 5.5) for all combinations of initial velocity 
and failure times of practical interest. Such an approach is impractical 
because of the number of combinations that would need to be considered. 
Rather, an approximate method of modifying the existing data (Fig. 5.4 


and 5. 5) to include the initial conditions must be found. 


Two approximations are made in this regard, and the errors 
induced in practical problems are evaluated. First, the failure time condi- 
tions is approximated by applying to the debris particle a negative impulse 
equal to the area under the dynamic pressure time curve from time t= 0 
to time t = t,. Mathematically then, the apparent initial impulse is given 


f 
in nondimensional form as, 


- 
2 2,2 
pa Cy, t “q < 
Gis gyoy- ——S 9 9 fe 47 (1-7 ae 
ig m 
° 

where 

a = apparent initial velocity 

T1(0) = actual initial velocity. 


when this is carried out the apparent initial velocity becomes, 


2 -2T 
ee e oe 1 1 
o. = G0) - oa - (tT, - Tr + sz ) - | Ps 


-2T 
The failure time Tr is always small (tT, <1) so that if e is expanded 
in a power series and ve is neglected in comparison to Te Then, 
92 
1 ~ 1 = 
t= b(0)- —> T, (5. 16) 


results. 
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The second approximation made was to assume that the particle 
motion resulting from the apparent initial velocity acting alone can be super- 
posed on the particle motion which is due to the blast wind, to give the re- 
sulting motion of the particle which is due to the combined effect of initial 
velocity and blast wind. By virtue of the general dependence of air friction 
on the square of the relative velocity, the two results are superposed by 
the square root of the sum of the squares. The motion which is due to the 
initial velocity can be determined from Eq. (5. 7) subject to the initial con- 


ditions 


6(0) 
G*(0) 


iT 
we 


The solution of Eq. (5.7) is then, 


ie fre 
1° tome 

rt! ok (5. 17) 
Oy = -k log = 


For the case of interest, 


In summary, the procedure for applying these results to a real 
problem is as follows: 
(1) Compute an from Eq. (5. 16) 


(2) Compute oY (1) and 6, (7) from Eq. (5.17) for the range 


of t's of interest. 


(3) Use Fig. 5.4 or 5.5 to determine o5(7) and 6 4(7) 


for zero initial conditions for the same values of T as above. 
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(4) At each given ( tT) compute the final results, 
P14)? + 04)? 


The validity of these results was tested by comparing the 


es 
! 


C' 


approximate results obtained by the above procedure to "exact" solutions 
obtained by numerically integrating Eq. (5.4) subject to the initial conditions. 
These comparisons are presented on Fig. 5.9. It can be noted that the re- 
sults are reasonably accurate in view of the over-all accuracy requirements 
of the debris prediction problem. Forvery high initial velocities (e.g., 
velocity displacement curves convex downward with a high initial peak) 

the approximate analysis is not very good. However, this problem is 


probably not very important when considering the nuclear environment. 
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CHAPTER SIX 
VULNERABILITY OF FIELD TROOPS TO TREE DEBRIS 


A supplementary study involving vulnerability of field troops to 
casualties from tree debris caused by a nuclear explosion in the proximity 
of a forest, instituted upon recognition of needs by the Office of the Chief 
of Engineers for estimating the hazards to engineer and field troops, was 
undertaken. Making certain simplifying assumptions (i.e., zero-strength 
tree limbs, plane blast wave loading, unobstructed trajectories, and that 
a hit upon personnel by a tree limb is a cagualty), tree limbs are followed 
in their trajectories from the time of shock loading to their impact with the 
ground. The safe distance may fall either inside or outside the forest and 


both cases are treated. 


Results of this analysis show that for the lower yields (1 KT, for 
example) a uniform horizontal translation of all branches is obtained since 
trajectories become vertical with an attendent small vertical drop, and the 
appearance of the area in front of the forest up to the ''safe distance'' would 
be similar to that of a forest floor after all of the branchwood were allowed 
to drop vertically. Basically similar behavior is observed for the higher 
yields where trajectories terminate before they become vertical (20 MT, 
for example), with the exception that the highest branches of the first few 
rows of trees pile up in a lower density than those following the closer-in 


trajectories. 


6.1 Previous Studies 


The results of previous studies on tree vulnerability which have 
influenced our assumptions concerning debris (Ref. 21, 22, and 23) are 


as follows: 


(1) Results of OPERATION UPSHOT-KNOTHOLE indicated 
that stands of 145 Ponderosa pines of heights 50 to 
75 ft offered no attenuation of peak overpressure or 


dynamic pressure. 
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(2) That low burst heights were found to cause more 
damage to trees than large burst heights when 
the peak dynamic pressure was the same. 

(3) OPERATION CASTLE indicated no pressure attenua- 
tion from the trees in natural tree stands. Damage 
predictions for two weapon yields compared favor- 
ably with the observed damage. 

(4) Damage to broadleaf stands is principally limb 
breakage and defoliation with occasional breakage 
of the main stem or uprooting. 

(5) The deflection and breakage of trees in the stand on 
UPSHOT -KNOTHOLE Shot 9 were approximately 
twice the values predicted on the basis of calculations 
of the first maximum deflection. By including the 
probability of breakage during the second maximum 
deflection under the negative phase, the predictions 
were brought into agreement with experimental values. 

(6) Trees are drag structures. The best parameter 
with which drag can be correlated is the dry weight 
of the crown. 

(7) In OPERATION SNAPPER, it was observed that when 
stem breakage occurred the stems broke at the 
tree base. 

(8) A considerable amount of data describing the 
mechanical and aerodynamic properties of tree stems 
and crowns is available. Almost none exist for the 


isolated branches. 


6. 2 Problem Approach 


Consider first the case shown in Fig. 6.1 where the troops are 


dispersed outside of the forest. 
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Figure 6.1 Relative Positions of Troops, Forest, and Ground Zero 


Our objective is to determine a conservative or upper bound ''safe distance'' 
rather than tackle the enormous task of finding the ''smallest safe distance''. 


To accomplish this goal, we made the following assumptions: 


(1) Ground burst is assumed for all attack conditions 
since it produces the most severe tree damage. 

(2) Trees are assumed not to attenuate static or dynamic 
pressures. 

(3) The trajectories of the flying branches are considered 
to be unobstructed by the other branches and trees. 

(4) If any debris strikes an individual, he is assumed to 
be incapacitated. 

(5) The total flight time" of a branch is considered to 
be equal to its free-fall time ina vacuum. This as- 
sumption turns out to be unimportant for yields no 


greater than 1 MT. 
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(6) Positive wind phase is taken equal to the positive 
shock wave phase. They are approximately equal. 
(7) The branches travel the same horizontal distance 


as an air particle during the positive wind phase. 


Assumption 7 is the most far-reaching and important on the list. 
Physically 1t corresponds to a branch which has zero strength, is complete- 
ly diffraction insensitive, and which has an infinite acceleration coefficient. 
The acceleration coefficient is the product of projected area and drag co- 
efficient divided by mass. Clearly, for tree branches this is a relatively 
large factor. Assuming the branches to be only drag sensitive is fairly 
good; assuming them to have zero strength is conservative but poor. The 


items tend to counteract one another. 


If only drag forces were operative, it is quite clear that a 
particle of air would be transported further than a solid object. The 
questions arise when the solid object has an initial velocity; for example, 
particles originating as crater throwout. Such particles are not of concern 
here; however, it is conceivable that sufficient impulse 1s delivered to a 
branch to both sever it from the tree stem and give it an initial velocity 
greater than the peak particle velocity. In such a case the branch or any 
other particle would experience a deceleration resulting frorn a drag force 
opposing its motion. The high acceleration coefficient of a branch would 
quickly bring its velocity into coincidence with the particle or wind velocity 
behind the shock front. 


During the negative phase, ihe reversed winds decelerate any 
airborne objects. Particles haveing sufficiently high acceleration 
coefficients have their forward motion reversed. In the spirit of conserva- 


tism, the negative phase is neglected. 


At low pressure levels (2.4 psi), the stems of trees remain 
standing and offer considerable interference to the flight of branches. 
However, at pressures of interest (between 5 and 30 psi) the stems are all 
broken and probably create no interference to flying debris. Because the 
shock wave moves at a higher velocity than the winds behind it, there isa 


strong possibility that debris from remote trees travels ahead of debris 
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from trees closer in to ground zero. 


The trajectories of the wind particles have been computed as de- 


scribed in Appendix H and the results are presented in Fig. 6.2. The fol- 


lowing remarks are based on, or are concerned with, this figure. 


(1) 


(2) 


(3) 


(4) 


(5) 
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Pressures over 30 psi are not included in trajectories 
since these pressures represent a greater hazard than 
the debris. 

As the wind velocity drops to zero the traje tories 
become vertical. The curves for the 20-MT devices 
were terminated at about 255 ft of vertical drop since 
trees of greater height are not of interest. 

The horizontal distance associated with the vertical 
tangent to the trajectories can be scaled approximately 
according to cube-root scaling; e.g., D/D, = wiw))!3 
where D 18 the distance and W the yield. 

In the cases where the trajectories become vertical 
with an attendant small vertical drop (e. g., the case 
of the 1-KT device), a uniform horizontal translation 
of all the branches was obtained. This situation is 
illustrated in the sketch shown in Fig. 6.3. The 
appearance of the area in front of forest labeled ''safe 
distance from forest'' would be similar to that of a 
forest floor after all of the tree branchwood were 
allowed to drop vertically. This situation is illustrated 
in Fig. 6.4 and 6.5 for forests exposed to low pressure 
blasts. 

In cases where the trajectories terminate before they 
become vertical (e.g., the case of the 20-MT device), 
basically the same behavior as in the previous case was 
found, with the exception of the highest branches 

on the first few rows of trees. These branches pile 

up with a lower density than those following the closer- 
in trajectories. This situation is depicted in Fig. 6.6 


where the final horizontal locations of branch centers 
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Figure 6.3 Trace of Branchwood for Low Yield Weapons 


for various branch elevations in a typical forest are 
indicated. Each cross represents the branchwood 
in a 5-ft vertical distance along the tree stem. 

The most cursory examination of typical trees 
indicates that the limbs in the top 15 ft would com- 
pletely cover the area surrounding the tree. Refer- 
ring to Fig. 6.6, we find that "complete kill" is ex- 
perienced for distances up to 1800 ft. Going from 
1800 ft to 1880 ft, the density of limbs diminishes 
and 100 percent kill is not expected in this region. 
The precise determination of the kill probability in 
this 80-ft area is not warranted and in such case the 
safe distance from forest is set equal to the maxi- 


mum trajectory (e.g., 1880 ft). 
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Figure 6.4 


Forest Stand After a Nuclear Explosion 


(2. 4 psi overpressure) 


Figure 6.5 


Forest Stand After a Nuclear Explosion 


(3.8 psi overpressure) 
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Figure 6.6 Horizontal Displacement of Tree Limbs 


We now turn our attention to the very straightforward task of 
determining the safe distance from ground zero in a forest of infinite extent. 
In Table 6.1 we have reproduced a schedule of damage criteria for forests 
from the ART 6. 24 of Reference 20. The damage level described in class D 


of this table represents the borderline condition for troop safety. 


Table 6.1 
DAMAGE CRITERIA FOR FORESTS 


Damage Nature of Damage Equivalent Hurricane 

Class Wind Velocity, 
(mph) 

AandB Up to 90 percent of trees blown down; 130-140 


remainder denuded of branches and 
leaves. (Area impassable to vehicles 
and very difficult on foot. ) 


Cc About 30 percent of trees blown down; 90-100 
remainder have some branches and 
leaves blown off. (Area passable to 
vehicles only after extensive clearing.) 


D Very few trees blown down; some leaves 60-80 
and branches blown off. (Area passable 
to vehicles. ) 
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Table 6.2 
SAFE DISTANCES TO PREVENT CASUALTIES FROM TREE DEBRIS 


Safe Distance Falls 


Yield Inside of Forest ” Safe Distance from Ground Zero 
of Sos for Various Overpresoures for Various Overpressures 
Weapon Safe Distance at the Front of the Forest, at the Front of the Forest, 
from Ground Zero (yd) ‘ (ya) 
[Spa] 10 pat [15 pel] 20 poll 30 pal |S pel] TO pal] 15 pol [0 pat [0 pal 
98 87 7 
100 
0.05 KT 376 yards 
0.1 KT 473 yards 
O.5 KT 809 yards 


KT 845 - 1021 yards 
{ 

Joes 

! 20 5,119 

40 5,137 

60 5,147 

80 5,153 

MT 100 8,448 - 10,208 yards 5,157 

120 5, 160 

160 5,163 

200 5, 163 

20 13,792 

40 13, 821 

60 “113, a42 

| 80 13, 858 

LomtT | 100 22,932 - 27,709 yarde | 13,871 

120 13, 882 

160 13,901 


L 200 13,916 
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The wind velocities 60-80 mph are associated with overpressures of 1.7 - 
2.2psi. The distance from ground zero for a surface burst at which these 


pressures are realized can be scaled from Fig. 3. 94a in Reference 20, 


6.3 Results 


Based on the approaches described in this chapter, we have 
computed the safe distances to prevent casualties from tree debris for 
various conditions and presented them in Table 6,2. Whenever the safe 
distance from the forest, which was computed from the maximum trajectory 
range, fell below the associated tree height, the tree height was used as 
the safe distance from the forest. At the pressure levels considered the 
trees will be blown over; hence, the stems and not the branchwood represents 


the more severe hazard under these conditions. 
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CHAPTER SEVEN 
VULNERABILITY OF FIELD TROOPS TO THROWOUT DEBRIS 
FROM CRATERING AND STREAM-BED CHARGES 


This study was undertaken to define a ''safe line" for positioning 
engineer or field troops in the proximity of very-low yield nuclear cratering 
and stream-bed charges, based on debris criteria. Since no reliable 
analytical method of predicting crater throw out debris was available, 
the problem involved locating and utilizing experimental data on debris dis- 
tribution under variations in the major controlling parameters - weapon 
yield, depth of burst, and soil characteristics. A number of sources were 
found to include data concerning crater throw out debris, some including 
variations in parameters (Ref. 2, 24, 25, 26, 27, 28, 29). Of these, two modes 
of measurement are used: the U.S. Geological Survey (Ref. 24) expresses 
debris distribution in terms of fragment sizes; the Suffield Experiment 
Station (Ref. 27) presents debris data in terms of a real density; and the 


Boeing Airplane Company (Ref. 25 and 26) presents data in both ways. 


The procedure described here for estimating the throwout 
environment about cratering charges is based primarily on the U.S. Geo- 
logical Survey reports of cratering tests in basalt in Area 16 at the Nevada 
test site (Ref. 24). Contours were plotted by USGS defining ground ranges 
for several average particle sizes. Average ranges for various fragment 
sizes were derived from these plotted contours and related to the yield 
and depth of burst. Nomographic methods are included whereby the average 
debris distance can be obtained for a wide range of weapon yields, depths 
of burst, and fragment size in basalt. A multiplication factor is introduced 
to convert the average ground range to maximum ground range, based on 
observations of the ray-like patterns noted in the USGS report on basalt 
craters. Data from the earlier Panama Canal series of tests were used 
to develop correction factors in converting the estimates for basalt to 
estimates for other soil media (Ref. 30). Likewise, debris measurements 
for cratering tests in marine muck, conducted as part of the Panama Canal 
series of experiments, were used to provide estimates of debris distance 


for streambed charges. 
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71 Method of Solution 


In general, the debris beyond the crater lip has been observed 
to conform to the following two expressions; 


1 
fod (7.1) 
x 1 
and 
Cc 
2 
Ds (7. 2) 
x2 
where 
p = areal density in weight of debris per unit area 
= fragment size 
x = distance from ground zero 
Cc) C, = constants 
Nyen, = exponents 


Data which verify Eq. 7.1 are available (Ref. 2, 25, 26, and 27). 
Data which verify Eq. 7. 2 are also available (Ref. 2 and 26). 


Equations 7.1 and 7. 2 comprise the basis for two different sets 


of relationships which can be used for determining safe line distances. 


Certain question regarding the use and validity of Eq. 7. 2 must 
remain. First, there is the greater tendency for fragments which travel 
greater distances to break into smaller fragments when they hit the ground. 
Therefore, even if all fragments are of equal size at takeoff, a decrease in 
particle size for increased distance would still be witnessed. Secondly, 
when air resistance is considered, it can be reasoned that, in general, 
large particles will travel further than small particles. Of course, there 
are other considerations such as the origin of different-sized particles 
relative to the point of burst, which lend reasonability to the observed 
distribution. It is not known which consideration is most important. In 
future tests, it would be desirable to consider the occurrence of impact 
breakage, since the criteria for safety involves the fragment characteristics 
before impact rather than after impact. 
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Let us now consider the criteria for the determination of a safe 
line ground range. The Armed Services Explosives Safety Board has used 
a 1-lb fragment as being capable of producing a fatal injury. If the safety 
criteria were established considering any injury less than fatal as safe, 
then if material density is known it is possible to use Eq. 7. 2 to find the 


safe line ground range. 


To apply the safety criteria to Eq. (7.1), consider the area 
covered by an average soldier laying flat on the ground. For a given areal 
density, the most severe situation is that in which all of the material landing 
in the area covered by the soldier were lumped into a single fragment. The 
distance at which the areal density is such that would require this fragment 


to be equal to the critical weight (e.g., 1 lb) is the safe line ground range. 


The one-pound fragment considered here is only an example. 
The actual size of fragment selected is a function of the probability of injury 
which is considered desirable. A much smaller fragment size might be 


selected if velocities were high. 


When a criterion for determining the equipment damage capa- 
bilities of fragments is developed, it will also be possible to apply Eq. (7. 1) 
and (7. 2) to equipment. 


Before any problems can be solved, we must find the functional 
relationship between the constants C and n, and the independent blast 
parameters. These parameters consist of yield, depth of burial, and soil 
characteristics. A fourth parameter, geological conditions, is also impor- 
tant (Ref. 24). It is very difficult, however, to evaluate geological factors, 
especially in military situations, and these will therefore not be considered 
as aparameter. Let us first attempt to develop relationships between the 
three blast parameters and the C and n constants for the distribution 


in fragment size. 


7.2 Fragment Size Distribution Method 


Ross B. Johnson (Ref. 24) has used maps to present data concern- 
ing throw out from a series of high explosion craters in basalt. Contour 


lines of equal average particle size for 1. 0-ft and 0. 5-ft diameter particles 
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are presented. (See Fig. 7.1, for example.) This information is used to 
obtain a first approximation of the relation between C, n, and the blast 
parameters. It must be emphasized at this point that the relations developed 


here are only first approximations. 


Data from eight charges of 1000 lb TNT and three charges of 
40, 000 lb TNT exploded at various depths are tabulated in Table 7.11. 
The areas within the contour lines of each map were measured and the radii 
of circles of equivalent areas calculated. Each radius so calculated is 
termed the "area mean distance for particles of given size". Accepting 


the size distribution law, 


erie (7. 3) 


where Eq. (7.3 is simply a modification of Eq, (7. 2), the area mean distance 
was plotted as a straight line on a log-log plot (Fig. 7.2). The slopes of the 
curves n and the coefficients C were found and tabulated in Table 7. 1. 
Table 7,1 expresses C in scaled form. The resulting data were then used 
in plotting Fig. 7.3 and 7.4. Figure 7.3 is a plot of the Distribution Law 
Exponent, n, vs scaled depth of burial, asw!/3. 7 and Fig. 7.4 is a 

plot of scaled distribution law coefficient, c/w! /3. .. vs scaled depth of 
burst, a/wi/3-5, 


A scaling factor of 3. 5 was found to produce the best fit. Its use 
agrees with data presented by R. B. Vaile of Stanford Research Institute 
(Ref. 31). Figure 7.5 shows the Vaile curves of crater radius vs yield for 
different materials. The curve for sandstone yields a scaling exponent of 
3.6. Since Johnson's data are for basalt it may be expected that they also 


have ascaling coefficient in the vicinity of 3.6, hence, 3. 5. 


Note that both curves of Fig. 7.3 and 7. 4 have the general appear- 
ance of an inverted parabola, similar to the curves of scaled crater radius 
and scaled crater volume versus scaled depth of burial found in many 
references on cratering. Notice on Fig. 7.4 that not all points fit the curve 
closely. At first glance one might expect an elipse. However, the results 
of the 1000-lb TNT shots are encouraging. Each data point for the 1000-1lb 


shots represents two shots whereas each point for the 40, 000-lb shots 
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Figure 7.1 


High Explosion Crater 13 (40,000 pounds), Buckboard Mesa, 
Nevada Test Site, Nye County, Nevada 
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Post-explosion topographic contours (by Am. Arial Surveys) 


Interval 10 feet. Datum is sea level 
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Figure 7.2 Area Mean Distance versus Fragment Size 
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Figure 7.5 Crater Radius versus Yield, Surface Shots, Various Soils 
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represent only one shot. Data for the 1000-lb shots show in general, what 
would be expected, i.e., a curve similar to an inverted parabola. Especially 
encouraging is the fact that the vertex of the parabola lies very close to the 
optimum scaled depth of burial for basalt. The data for the 40,000-lb shots 


show a great deal of scatter. 


Figure 7.6 is used to determine the equations for the curves of 
Fig. 7.3 and 7.4. Offsets from axes placed through the vertical were 
plotted in Fig. 7.6 to obtain coefficients and exponents. The resulting 


equations are: 


n= 1,28- 0.25 (A, - mea (7. 4) 

C= wr? E 34 - 0.52 (A, - 2. 2373] (7. 5) 
where 

n = distribution law exponent 

C = distribution law coefficient 

W = charge weight in lb of TNT 

de - aswi/3- 5 


d = depth of -turial 


Equations (7.3), (7.4), and (7. 5) are the basis of the charts of Fig. 7.7 and 
7.8. These charts can be used to predict safe line ground range. Let us 
reiterate that this represents only a first approximation. Instructions for 


the use of the charts follows. Two examples are presented on the charts. 


7.2.1 Nomographic Calculation of Safe Distances 
Based on Size Distribution 
Instructions for Use of Charts 
Find n, and C as follows: 


(1) Select yield of device in Fig. 7.7 on horizontal axis. 
(2) Draw a vertical straight line upward to the diagonal 


line representing the desired depth of burst, d, (ft). 
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Values of Distribution Law Coefficient and Exponent 
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Figure 7 € Distribution Law Exponent and Constant 


versus Scaled Depth of Burial 
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(3) Draw a horizontal line from the resulting intersection 
with the diagonal to both the n and C curves. 

(4) Draw a line vertically upward from intersection on n 
curve to the top scale and find the value of the distri- 
bution law exponent n. 

(5) Draw a line vertically from the intersection with 
the C curve to the line A-A., 

(6) From the intersection with A-A, draw a horizontal 
line to intersection with the diagonal (parallel to 

A-A) line of constant device yield. 

(7) Draw a vertical line from the resulting intersection 

upward to the C scale and find the value of the distri- 


bution law coefficient, C., 


Use values of n and C found above as input for Fig. 7. 8 and 


find maximum average distance of encounter with particle of diameter D. 


(1) Select diameter of particle on the lower horizontal axis. 

(2) Draw a vertical line upward to the diagonal represent- 
ing the value of n found previously. 

(3) Draw a horizontal line from the resulting intersection 
to the diagonal line of constant value of the C found 
previously. 

(4) Draw a line vertically to the upper horizontal scale to 
find the maximum radial distance from ground zero at 


which the particle of diameter D may be expected. 


Note now that distances found by this chart are average ground ranges at 
which a given average particle size will be found. In reference to Fig. 7.1 
it will be recognized that a large degree of variability can be attributed to 
the geological characteristics of the soil, such as discontinuities, non- 
homogeneity and anisotropy. 

To estimate a safe line, it is necessary to take these factors into 
consideration. This can be done by considering the maximum variation 
about the average. Table 7. 2 tabulates the variations by showing the calcu- 


lation of the ratio of maximum contour line distance and average contour 
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line distance. The resulting ratios were plotted on semilog paper versus 
yield in lb of TNT in Fig. 7.9. 


Note the general decrease in x Ix ratio for increasing 
max’ “avg 

yield. This is exactly what would be expected for the following reason. 
The soil will appear to be very nonhomcgeneous, discontinuous, and 
anisotropic to a small blast. As the blast yield increases the appearance 
of the soil will become homogeneous, continuous, and isotropic. The 
radomeness in throw out will therefore become less as blast yield is 
increased and consequently we can expect the asymmetry ratio-vs-yield 


curve will be asymptotic to the line 


a 


max _ 


Data from DANNY BOY have been included in Fig. 7.9. However, 
compatibility with the other data is questionable since these values of 
asymmetry are based, at least in part, on elevation contours instead of 
particle size contours. The term asymmetry is used since it describes the 
throw out pattern. Figure 7.10 shows the DANNY BOY contours. 


Figure 7.9 can be used to find an asymmetry factor. This factor 
should be multiplied by the average contour line distance found from Fig. 7. 8. 
Since the results are still only a first approximation it is advisable to apply 
a safety factor of at least 1.5 to the results. Additional charts can be de- 
veloped to take the asymmetry and safety factors into account automatically. 


The final result will be a safe line ground range for blast set off in basalt. 


If it is desired to scale from basalt ot the other materials, 
Fig. 7.11 will provide factors which can be applied to the final result obtained, 
However, in some cases such as Cucaracha and Culebra, and inarine muck, 
this scaling may be questionable since it is so difficult to consider fragments 
when speaking of these soils. The scaling may possible be good only for 
the materials capable of fragmentation. The curves of Fig. 7.11 were 
developed from the debris curves of Fig. 7.12 taken from the Panama Canal 
Studies (Ref. 29). 
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Cucaracha and Culebra is 
a very weak Bentonite 
material, moist, with sand 
and gravel mixed in 


1/3) 


Figure 7.11 Debris Diameter Ratio versus Scaled Depth of Burial 


SECRET 


201 


SECRET 


Scaled Debris Diameter, ( {t/lbs! 13 ) 


Scaled Depth of Charge, ( ft/lbs! 13 ) 


Figure 7 12 Debris Diamete- for HE Charges in Various Soils 
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7.2.2 Explosive Equivalence 


The data presented in Fig. 7.7 and 7.8 are for TNT charges. 
If nuclear charges are to be used, it is necessary to find an equivalence 
factor which would provide a means of finding the TNT equivalent of the 
nuclear device. The state-of-the-art concerning equivalence is uncertain 
at present. One can only hypothesize on the basis of the limited data 


that exist- 


The first problem concerning equivalence is one of semantics. 
Depending upon the terminology used, almost any factor ranging from 2 
percent to 120 percent is acceptable. First, let us here define equivalence. 
For this study, a TNT and a nuclear device would be equivalent when each 
produces the same amount of throw out debris distributed at the same 


distance. 


One approach to equivalence is to hypothesize on the basis of the 
energy available in a nuclear blast. Effects of Nuclear Weapons (Ref. 20) 
states that energy of a nuclear blast is divided as follows*: 


Blast and Shock - 50 percent 
Thermal Radiation - 35 percent 
Residual Nuclear Radiation - 10 percent 
Initial Nuclear Radiation - 5 percent 


On this basis it can be expected that at the ground surface a nuclear blast 

will be at least 50 percent efficient in terms of debris-producing capabilities. 
If the blast is set off underground a portion of the 35 percent thermal radiation 
will be transformed into mechanical effects as a result of the vaporization 

of material in the immediate vicinity of the device. Therefore, for buried 
bursts efficiencies between 50 percent and 85 percent are expected. The 
above does not take into account variation due to partition of energy between 
the air and soil - this could, in fact, result in a substantially lower 

efficiency. For our purposes the more conservative estimate has been 


selected. 


* Taken in atmosphere. 


SECRET 


203 


SECRET 


For bursts near the surface, venting will occur and therefore 
a portion of the thermal radiation will still escape. On this basis one would 
expect a curve of efficiency versus scaled depth of burst te be similar to 


that of Fig. 7.13 


1.00 
0 85 


Efficiency 
° 
131) 
° 


ajwi/3 


Figure 7.13 Percent of Efficiency vs Scaled Depth of Burst 


7.2.3 Comparison with DANNY BOY Results 


Maximum distances for fragments of various sizes have been 
computed by the methods of this chapter, using DANNY EOY explosion 
parameters as inputs. This was an 0.430-KT device buried at 110 ft in 
basalt. Predicted debris distances are plotted in Fig. 7.14 for various 


equivalence factors according to the methods described here. 
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100% Equivalence 
85% Equivalence 
1100 


50% Equivalence 


eweweee 25% Equivalence 


Note: Assymetry factor of 1, 2 was 
used, no safety factor included. 
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Figure 7.14 Estimated Maximum Fragment Distances for DANNY BOY Event 
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7.3 Throwout Debris from Stream-Bed Charges 


Since energy partitioning for an explosion at the interface 
between two media is inversely proportional to the ratio of the densities 
of the media, the bottom burst of a stream-bed charge would be expected 
to impart more blast energy to the soil than a comparable surface burst 
at an air-soil interface. Thus, the stream-bed charge would dislodge 
more material for contribution to the debris, which we shall here regard 
only as the above-water ejecta. Initial velocities of the ejecta from 
stream-bed charges are unknown, but on the basis of energy partitioning 
we can concede, conservatively, that they may exceed initial velocities 
of ejecta from the comparable surface burst on land. Ejecta from the 
bottom burst may follow three paths, assuming the initial gas bubble does 


not break the surface. 


(1) Upward (nearly vertical) within the rising gas 
bubble into the atmosphere at relatively high 
velocities. 

(2) Upward (nearly vertical) through water, with 
considerable retardation. 

(3) Through the water at various angles to the 


water surface, with considerable retardation. 


Bottom charges in deep stream beds would not be expected to 
throw debris to greater distances than comparable surface bursts on land, 
for only the debris rising through the gas bubble would be expected to reach 
the atmosphere at velocities comparable to those of the ejecta from surface 


bursts on land. 


Although little consideration has been given crater throwout 
debris from underwater explosions or bottom bursts, either nuclear or high 


explosive, the following observation supports the above conclusions. 


(1) The displaced bottom material did not produce air- 
borne debris considered of consequence in the ''Baker" 
test of OPERATION CROSSROADS. This event involved 
an explosion of a 20-KT device at 180-ft depth in the Bikini 
Lagoon, which is considred relatively shallow for that 
yield (Ref. 30). 
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(2) Photographic observation of a series of tests in- 
volving underwater and bottom bursts of 300-1lb 
HE charges showed plumes with decidedly vertical 
trajectories, (Ref. 31). This behavior is shown in 
Fig. 7.15 (from. Ref. 31) and 7.16 (from Ref. 20). 

(3) Full-scale nuclear test experience has shown that 
if the depth of the underwater burst is not too great, 
the bubble remains intact until it rises to the water 
surface. At this point debris is expelled into the 
steam and fission gases (Ref. 29). Such debris 
would have predominantly vertical trajectories as 


shown in Fig. 7. 16. 


The problem of far-flung debris from underwater bursts is 
probably a hazard to field troops in cases of shallow streams only -- where 
the initial gas bubble breaks the surface and provides ejecta entering the 
atmosphere at all angles and with high initial velocities. Under these 
conditions the stream-bed charge placed either on or somewhat beneath 
the bottom -- on the basis of energy partitioning -- can conceivably produce 
ejecta which enters the atmosphere at initial velocities greater than that 


from the comparable surface or shallow-burried burst on land. 


An analytic solution for initial velocities of the ejecta for either 
of the cases cited is beyond the scope of this portion of this investigation, 
and in fact, may be well beyond the scope of current knowledge. Experi- 
mentation to date can only contribute to rudimentary estimates of debris 


behavior from stream-bed charges. 


For explosions where the gas bubble does not initially break 
the surface or approach the surface, we can probably neglect maximum debris 
distance in determining a safe line for personnel. In these cases radioactive 
spray from the plume itself, the condensation cloud and the base surge may 


extend farther than the ejected bottom material. 
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Figure 7.16 
Plume From Shallow Underwater Nuclear Explosions 
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Among all the experimental work reviewed on cratering and debris, 
the Panama Canal tests of cratering charges in marine muck most closely 
resemble the configuration of the shallow stream-bed charge -- at least 
they assuredly provide a situation where the ejecta is not retarded in its 
trajectory by passing through water. Energy partitioning may not be 
exactly comparable to the stream-bed charge, for the immediate escape 
of the gas bubble would also cause more energy partitioning away from 


the bottom material. 


In the absence of analytical methods for estimating debris behavior 
from bottom bursts, recourse is made to the available experimental data on 
crater tests in marine muck. Debris limits for stream-bed charges can be 
computed by estimating debris distances for basalt as prescribed earlier 
for cratering charges, and then applying a correction to account for the 
ratio of debris limits for basalt and bottom material. This correction 


factor, 


Debris diameter for bottom material 


“Debris diameter for basalt 
can be taken as the debris-diameter-ratio for marine muck in Fig. 7.11. 
In using Fig. 7,11 in this manner the scaled depth of burial of the charge 
should be taken relative to the stream bed. Thus for bottom bursts, scaled 


depth of burial will be zero. 
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APPENDIX A 
MAXIMUM DEBRIS DISTANCE AND EXPLOSION PARAMETERS 
FOR SELECTED EXPLOSIONS 
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MAXIMUM DEBRIS DISTANCE AND EXp 
FOR SELECTED EXP 


re o 2S 

oo E 8 gE 

re 62 2a ‘a3 

w Zo on ¢2 ec og sa 

: g | Rg x5 se6 | 222 

wa > als EO ag2 aac 
Indiana Ordnance Wks Smokeless Powder | Building 81.59 0.58 
Beira, Portugal Black Powder Mag. 82.50 0.58 
Highland Station, Calif. Dynamite Unbarricaded Building 103.6 0.79 
Haskell, N.J. Smokeless Powder | Unbarricaded Building 117.0 0.58 
Rio de Janeiro, Brazil Dynamic Unbarricaded Building 124.0 0.79 
Arco, Idaho 50/50 Amatol Barricaded Igloo 125.0 0.87 
Savanna Ordnance Depot Tetrytol Earth Covered Igloo 147.0 1,20 
Kobe, Japan Dynamite Barge 150.8 0.79 
Manila, P. 1. Dynamite, etc. Unbarricaded Mag. 50x150 173.1 0.79 
Indiana Gun Powder Frame Bldg. and RR Cars} 177.6 0.58 

Rifle Powder 
Charleston Rifle Powder Magazine 179.2 0.58 
Tessenderloo, Belgium Ammon. Nitrate Unbarricaded Building 193.0 0.42 
ae eee TNT Freight Cars 200.0 1,00 
Mindi Magazine, Canal Zone Dynamite Unbarricaded Mag. 225.5 0.79 
Sonemachi, Japan HE Unbarricaded Dump 270,0 1,00 
Baltimore, Md. Dynamite Steamer 300.0 0.79 
Acisate (Varese), Italy Ammo. Underground Bunkers 350.0 1.00 
Guadalcanal Torpex Steamer 400.0 1,25 
Bari, Italy Comp. B Steamer 544.0 1,10 
Hastings, Neb. Torpex, TNT and | Barricaded Bldg.500 x 25| 550.0 1,17 
DB Powder 
Pleasant Prairie, Wis. Black Powder Unbarricaded Mag. 587.5 0. 58 
Bombay, India Steamer 400.0 1,00 
Lake Denmark, N.J. Unbarricaded Mag. 800.0 1,00 
Bucharist, Rumania Building 1,000.0 1,00 
Mt. Hood, Pacific Theatre Steamer 1,000.0 1.00 
Brest, France Ammonium Nitrate | Steamer 730.0 0.42 | 

Texas City, Texas Ammonium Nitrate | Steamer 2, 280.0 0.42 
Port Chicago, Calif. Steamer 2,136.0 1,00 
Halifax, Nova Scotia Steamer 2,600.0 1,00 
Burton-on-Trent (Fauld) Misc. Bombs Barricade 1.00 
Oppau, Germany Ammonium Nitrate | Open Pile 0.42 
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MUM DEBRIS DISTANCE AND EXPLOSION PARAMETERS 
FOR SELECTED EXPLOSIONS 
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a | 3 E e 
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: S§ o22|au86 | 4 angles} so (as 
E >> Bei) 8§ssz )  & Ee | 481521 08 e gs 
oe na ad tt oak = Bee fae fn fe a mo. 8 
vin Eoa ORAL moe [m* a sha _|m' se | ss iB ao § 
& 8 125 Fen | Soe8 IL 8 $222 asi eeiFe [222 
EO ans daclasas [eS SAA~|BsA| 2A] ES Fac 
Building 8 0.5 43.75 3.62 900 249 15 11, 268 3,119 
Mag. 82.50 0.5 47.85 3, 62 2, 650 733 8 12, 385 3, 416 
Unbarricaded Building 0.7 81.80 4,34 7,920 1,828 8 15, 052 3, 472 
Unbarricaded Building 67.90 4.08 2, 250 552 15 12,991 3, 190 
Unbarricaded Building 98.00 4.60 3, 000 652 14, 976 3, 254 
Barricaded Igloo 108.8 4,77 3,950 828 15, 893 3, 335 
Earth Covered Igloo 176.3 5.60 6,000 [1,071 19, 023 3, 398 
Barge 119.0 4.91 | 17,920 | 3,650 15, 383 3,131 
Unbarricaded Mag. 50x150 136.9 5.15 1, 800 350 15,510 3,015 
Frame Bldg. and RR Cars 103.0 4.69 7,920 1,691 15, 268 3, 262 
Magazine 104.0 4.70 2,400 511 16, 241 3, 492 
Unbarricaded Building 81.05 4.32 5,280 |} 1,220 15 13,921 3, 221 
Freight Cars 200.0 5.85 5, 280 903 5 21,193 3, 630 
Unbarricaded Mag. 178, 2 5.62 7,920 | 1,410 8 19, 879 3, 539 
Unbarricaded Dump 270.0 6.46 2, 400 372 19, 650 3, 046 
Steamer 237.0 6.19 | 15,840 | 2, 560 20, 118 3, 257 
Underground Bunkers 75.0 4,21 1,500 356 24 12, 189 2, 894 
Steamer 500.0 7.93 9,000 | 1,135 20 26, 663 3, 366 
Steamer 599.0 8.42 4,500 535 20 28,498 3, 389 
Barricaded Bldg. 500 x 25 643.5 8.63 7, 300 846 30, 533 3, 544 
Unbarricaded Mag. 341.0 6.99 13, 200 1,890 23, 872 3, 424 
Steamer 400.0 7.36 3,900 530 20 24, 532 3, 336 
Unbarricaded Mag. 800.0 9.28 5, 280 569 15 32, 495 3, 508 
Building 1,000.0 10. 00 5, 280 528 15 35, 185 3,527 
Steamer 1,000.0 10, 00 6, 600 660 34, 395 3, 447 
Steamer 306.5 6.74 5, 280 784 22, 192 3, 297 
Steamer 958.0 9.85 11, 500 1,168 20 33, 857 3, 443 
Steamer , 136.0 12.88 13, 000 1,010 20 45,193 3, 518 
Steamer 2, 600.0 18, 480 20 50, 528 
Barricade , 670,0 13, 86 4, 346 314 20 48,924 3, 536 
Open Pile 1, 890.0 12, 34 4,920 399 35, 737 2, 898 
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Site 


Indiana Ordnance Wks 
Beira, Portugal 
Highland Station, Calif. 
Haskell, N.J. 

Rio de Janeiro, Brazil 


Arco, Idaho 

Savanna Ordnance Depot 
Kobe, Japan 

Manila, P.1. 

Indiana 


Charleston 
Tessenderloo, Belgium 
Black Tom Island, 

N. Y. Harbor 
Mindi Magazine, Canal Zone 
Sonemachi, Japan 


Baltimore, Md. 
Acisate (Varese), Italy 
Guadalcanal 

Bari, Italy 

Hastings, Neb. 


Pleasant Prairie, Wis. 
Bombay, India 

Lake Denmark, N.J. 
Bucharist, Rumania 

Mt. Hood, Pacific Theatre 


Brest, France 

Texas City, Texas 

Port Chicago, Calif. 
Halifax, Nova Scotia 
Burton-on-Trent (Fauld) 
Oppau, Germany 
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MAXIMUM DEBRIS DISTANCE AND EXPLOSION PAF 
FOR SELECTED EXPI.Qs) 


onfinement 
plosive 
quivalence 


Year 
Explosive 
Material 

T 

C 

Quantity of 


Smokeless Powder 
Black Powder 


Building 
Mag. 


O° | Ex 


wu | E 


Dynamite Unbarricaded Building 0.7! 
Smokeless Powder | Unbarricaded Building 0. 5! 
Dynamic Unbarricaded Building 0.7! 
50/50 Amatol Barricaded lgloo 0.8 
Tetrytol Earth Covered Igloo 1. i 
Dynamite Barge 0.7! 
Dynamite, etc. Unbarricaded Mag. 50x150 0.7' 
Gun Powder Frame Bldg. and RR Cars 0.5 
Rifle Powder 

Rifle Powder Magazine 0.5. 
Ammon. Nitrate Unbarricaded Building 0.4 
TNT Freight Cars 1,0: 
Dynamite Unbarricaded Mag. 0.7 
HE Unbarricaded Dump 1.0 


Steamer 
Underground Bunkers 
Steamer 

Steamer 

Barricaded Bldg. 500 x 25 


Dynamite 
Ammo, 
Torpex 


a) 
oe e @e e 
—eN ON 


Comp. B 
Torpex, TNT and 
DB Powder 
Black Powder 


oooow 


Unbarricaded Mag. 
Steamer 

Unbarricaded Mag. 
Building 
Steamer 


See Oo 
. e 


Steamer 
Steamer 
Steamer 
Steamer 

Barricade 
Open Pile 


Ammonium Nitrate 
Ammonium Nitrate 


Misc. Bombs 
Ammonium Nitrate 


0.4 
0.4 
1.0 
1,0 
1.0 
0.4 


MUM DEBRIS DISTANCE AND EXPLOSION PARAMETERS 
FOR SELECTED EXPLOSIONS 


» 0 —_— 
c V2 ea 
7) o; 2 2 
E 5 os 33 vb E ” 
Ww YM ont ot _ ad 13) 
Og aan | aoepok £8 G 
He O>H] >Ho Ba 
OH t= cm OD om 
a8 Sen | Bo oS ages 
al yaad ‘ast 
eG Age) seas SAA 
Building 0.5 
Mag. 0.58 
Unbarricaded Building 0.79 
Unbarricaded Building 0.58 


Unbarricaded Building 


Barricaded lgloo 

Earth Covered Igloo 
Barge 

Unbarricaded Mag.50x150 
Frame Bldg. and RR Cars 


Magazine 
Unbarricaded Building 


Freight Cars 


Unbarricaded Mag. 
Unbarricaded Duinp 


Steamer 

Underground Bunkers 
Steamer 

Steamer 

Barricaded Bldg.500 x 25 


Unbarricaded Mag. 
Steamer 
Unbarricaded Mag. 
Building 
Steamer 


Steamer 
Steamer 
Steamer 
Steamer 
Barricade 
Open Pile 
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15 11, 268 

8 12, 385 

8 15,052 

15 12,991 

15 14,976 


15, 893 
19,023 
15, 383 
15, 510 
15, 268 


oe Oe 
wm o Www 


16, 241 
13,921 


21,193 


19, 879 
19, 650 


~ 
own moe 


20, 118 
12, 189 
26, 663 
28, 498 
30, 533 


23, 872 


20 24, 532 
15 32, 495 
15 35, 185 


34, 395 


22,192 


20 33, 857 
20 45, 193 
20 50, 528 
20 48,924 


35, 737 
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APPENDIX B 
COMPUTER PROGRAM FOR REGRESSiON STUDY 
OF HE EXPLOSIONS 
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COMPUTER PROGRAM FOR REGRESSION STUDY OF HE EXPLOSIVE 
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APPENDIX C 
MAXIMUM DEBRIS DISTANCE AND CRATER DIMENSIONS 
FOR SELECTED EXPLOSIONS 
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Site 


Sunflower Ordnance Wks 
Perranporth, England 
Umbogintwini, Natal 
Sherborne, England 
Lower Hope Point, England 


Earle, N.J. 

Cabot, Pa. 

Alconbury, England 
Uplee's Marshes, England 
Selma, N. C. 


Barksdale, Wis. 
Honshu, Japan 

Ikego Ammunition Depot 
Winsted, Conn. 

Mt. Carmel, Pa. 


Badger Ordnance Wks., Wis. 
Ft. Belvoir, Va. 

Kurihama Naval Base, Japan 
Nebraska 

Cleveland, Ohio 


N.A.D., Oahu, Hawaii 
Railhead, Ordona, Italy 
McAlester, Okla. 
Marugama, Shikoku, Japan 
Sunflower Ordnance Wks 


Highland Station, Calif. 
Gibbstown, N.J. 
Chicago, Dl. 

Erith, England 

Kuba, Okinawa 


Manila, P.I. 
Antwerp, Belgium 
Reddick, Il. 
Black Tom Island, 
N.Y. Harbor 
Okinawa 
Johannesburg, So. Africa 
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M 
a 
CY 

‘al 


plosive 
aterial 


as 


Powder 
Gelatin 
NG 
TNT 
NG 


Torpex 

Black Powder 
Bombs 

NG 

Tetryl 


NG 

Bombs 

Ammo 

Black Powder 
Black Powder 


NG 

AN Cratering Explosive 
Ammo 

TNT 

Dynamite Black Powder 


Torpex 

TNT 

Black Powder 
Picric Acid 
NG 


Dynamite 

NG 

Black Powder 
Black Powder 
Dynamite 


Dynamite 
Black Powder 
Dynamite 


TNT - Picric Acid 


Dynamite 
Gelatin 


Quantity of 
Explosives, 
(tons) 


0. 250 
0. 266 
1.050 
0.250 
0. 500 


4.896 
3.000 
1.250 
0.50 

4.000 


4.835 
0.050 
74.00 
1. 250 
4.000 


3. 750 

0. 160 
32. 50 

8.085 
13, 375 


12.000 
6. 800 

22.47 
3. 750 
3.700 


103.6 
1.500 
81.00 

41,75 
1. 250 


173.0 
6.000 
15.00 


200.0 


48.01 
55.35 


1, 00) 
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quivalence 


(TNT 
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0.79 
1.06 


0.79 
1,35 


#%5 F 
Lore! ny 
ae au ue 
a 0 
soa 8 SES 
oY ¥ 6 ss oe 
aEaS eye) 
0.145 0.525 12.0 
0. 360 0.711 17.0 
1,418 1.122 18.0 
0. 250 0. 630 15,0 
0.675 0.876 16.0 
6.120 1,827 20.0 
1.740 1,201 25.0 
1,250 1.077 22.5 
0.675 0. 876 30.0 
4.560 1. 657 14,3 
6. 520 1. 867 40.0 
0.050 0. 368 30.0 
74.00 4, 200 30.0 
0.725 0. 898 30.0 
2.320 1, 322 40.0 
5.060 1.715 47.0 
0. 160 0. 543 39.5 
32.50 3. 188 40.0 
8.085 2.010 51.5 
9.235 2.097 60.0 
15.00 2.464 60.0 
6. 800 1, 894 45.0 
13.036 2. 372 50.0 
4.163 1, 607 75.0 
4.995 1,708 55.5 
8. 180 2.015 50.0 
2.025 1, 264 47.8 
46.98 3. 606 95.0 
24,22 2. 892 75.0 
0.988 | 0.995 75.0 
136.7 5. 150 120.0 
3,480 1,513 150.0 
11,85 2.279 120.0 
212.0 5.958 275.0 
37.92 3. 356 130.0 
74,72 4. 206 182.5 


}RIS DISTANCE AND CRATER DIMENSIONS FOR SELECTED EXPLOSIONS 
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1,590 
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3, 534 
4,245 
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53, 100 
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106, 000 
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265, 300 
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Crater Depth, 
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10, 980 
10, 615 


10, 700 
28, 400 
14, 730 
19,050 
55, 500 
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70, 000 
49, 500 


42,400 


79, 100 
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Distance, (ft) 
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5, 280 


1, 600 
9, 000 


SECRET 


Debris Disapces : 
Stnt 


3 
wi/ -Scaled 
(ft/tons 


SECRET 


APPENDIX D 
FRAGMENT-SIZE DISTRIBUTION AND DISPERSION 
OF FRAGMENTS 
FROM PANTEX ORDNANCE PLANT EVENT 
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FRAGMENT-SIZE DISTRIBUTION FOR REINFORCED CONCRETE STRUCTURE 
(Pantex Ordnance Plant Event) 


Cumulative 


Fragment Number Cumulative | Cumulative Weight Percent Cumulative 
Size of Number Percent of of Weight of Percent 
(iv) Pcacrmenta of of Total Fragments, | Total Fragments, of Total 
8 Fragments | Fragments (1b) Weight (ib Weight 
1/16 56 56 0. 18 3.5 0.0041 3.5 0, 004 
1/8 906 962 3.10 113,25 0.133 1 6.75 0.137 
1/4 18, 138 19, 100 61.47 4,534.5 5.324 4,651.25 5.461 
1/3 12 19,112 61.51 4.0 0.005 4,655. 25 5.466 
1/2 40 19,152 61. 64 20.0 0.023 4, 675, 25 5.490 
3/4 5,858 25,010 80.49 4,393, 5 5.159 9,068.75 10. 64 
1 32 24,042 80.59 32,0 0.038 9, 100.75 10, 69 
1-1/4 1 25,043 80. 60 1.25 0.001 9,102.0 10. 69 
1-1/2 10 25,053 80. 63 15.0 0.018 9,117.0 10.71 
2 2,502 27,555 88. 68 5,004.0 5.876 14,121.0 16.58 
2-1/2 8 27,563 88.71 20.0 0.023 14, 141.0 16, 60 
3 2, 680 30, 243 97.33 8,040.0 9.440 22, 181.0 26.04 
3-1/3 1 30, 244 97,34 3. 33 0.004 22, 184.0 26.05 
3-1/2 14 30, 258 97. 38 49.0 0.058 22, 233.0 26.11 
4 12 30, 270 97.42 48.0 0.056 22, 281.0 26.16 
5 2 30,272 97.43 t0.0 0.012 22,291.0 26,17 
6 13 30, 285 97.47 78.0 0.092 22, 369.0 26, 27 
7 1 30, 286 97.47 7.0 0.008 22, 376.0 26. 27 
8 3 30, 289 97.48 24.0 0.028 22. 400.0 26, 30 
9 1 30, 290 97.48 9.0 0.011 22, 409.0 26, 31 
10 7 30, 297 97.51 70.0 0.082 22, 479.0 26. 39 
12 106 30, 403 97.85 1,272.0 1.493 23, 751.0 27.89 
14 13 30,416 97.89 182.0 0.214 23, 933.0 28.10 
15 2 30, 418 97.90 30.0 0.035 23, 963.0 28.14 
18 9 30,427 97.92 162.0 0.190 24,125,0 28, 33 
20 1 30,428 97.93 20.0 0.023 24, 145.0 28,35 
40 4 30,432 "97.94 160.0 0. 188 24, 305.0 28,54 
50 1 30,433 97.94 50.0 0.059 24, 355.0 28. 60 
60 240 30, 673 98.72 14, 400.0 16.908 38. 755.0 45.51 
70 17 30, 690 98.77 1,190.0 1,397 39, 945.0 46.90 
75 41 30.731 98.90 3,075.0 3.610 43,020.0 50.51 
80 1 30,732 98.91 80.0 0.094 43, 100.0 50. 61 
90 248 30, 980 99.70 | 22, 320.0 26. 208 65,420.0 76.82 
95 2 30, 982 99.71 190.0 0,223 65, 610.0 77,04 
100 6 30, 988 99.73 600.0 0.705 66, 210.0 77,74 
150 1 30, 989 99.73 150.0 0.176 66, 360.0 77.92 
180 81 31,070 99.994] 14, 580.0 17,120; 80, 940.0 95.04 
225 1 31,071 99.997 225.0 0. 264 81, 165.0 95. 30 
4000 1 31,072 100.000 4,000.0 4.697 85, 165 100. 00 
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APPENDIX E 
FRAGMENTATION DATA 
ON EXPLODED DRY SANDSTONE BLOCKS 
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APPENDIX F 
CONCRETE FRAGMENT WEIGHTS AND DIMENSIONS 
FOR 1/24-SCALE SHIELDED REACTOR MODELS 
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CONCRETE FRAGMENT WEIGHTS AND DIMENSIONS 
FOR 1/24-SCALE SHIELDED REACTOR MODELS 


Event Simulated, 730 1b TNT, 1 msec, w/Pyrocore 
Vessel 3/4 Full of Water 


| Individual a Dist. of Fragments 
P 


Model No. 3 


iece Width, from Original Model R ks 
Number {| ( ) (in.) Position, ff >> 10 ft emarks 
12 Ring. 12-in. 
diam. with 
4-1/2-in. 
diam. hole 
43 ll} 2 1] 2 7 9 4 
1 8) 3 5 8 4 
4 8] 1 16 4 6 6 4 
2 5| 7 6 7 4 
6 5| 7 6 4 4 
5 5] 4 4 5 4 
8 5/1 21 3 6 5 4 
30 4165 4 5 6 4 4 
16 3] 11 4 5 4 
3 3] 9 4-3/4 4 4 
15 31 4 6 6 4 
32 3) 1 13 9 5-1/4 4 4 
2 2] 15 4 5 4 
9 214 4 4 4 
7 2133 4-1/2 4 4 
M1 2] 8 2-1/2 5-1/2 4 
37 2] 5 5 3 4 
24 2] 4 5] 11 4-1/2 3-1/2 4 
23 1} 15 4 3-1/2 3 
21 1] 14 5 2-1/2 4 14] 5 
17 1] 13 4 3 
10 1] 12 5 2-1/2 4 
20 1] 9 5 3 2-1/2 
22 1} 9 2-1/2 2-1/2 2 
27 1] 9 2 4 4 
28 1] 7 4 2-1/2 4 
31 1{ 7 4-3/4 2-1/2 4 
13 1] 5 5-3/4 2 4 
14 1] 5 2-1/2 4 4 
29 1 5 2 2-1/2 4 
19 1| 4 3 3 3 
33 1] 4 4-3/4 2-1/2 4 
18 ] 3 22 9 4 2-1/2 3 
40 15 4 2 2-3/4 
26 14 2 3-1/2 2-1/2 
36 13 4-1/2 1-3/4 3 
41 12 4-1/2 1-1/2 2 
25 1 2-1/2 1-1/2 4 
34 10 | 3/4 2 2-1/2 
38 10 4 1-1/2 3 
35 9 3-1/2 2 3 
39 7 4-1/2 1 2-3/4 
42 7 1-1/2 4 2 
45 6 1-1/2 2-1/4 1-1/2 18 
46 1/4) 3/4 ] 1/4 19 | 6 
4 1-3/4 | 1 4 Approximately 


| 50 small 
chunks of 
concrete 
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CONCRETE FRAGMENT WEIGHTS AND DiMENS:ONS 
Model No. 5 FOR 1/24-SCALE SHIELDED REACTOR MODELS 


Event Simulated, 510 1b TNT, 1 msec, w/Pyrocore 
Vessel Full of Water 


Individual , 
Weight, Madi Remarks 
_Weigh ( ) din.) 


hp bh hw ed ALL HH 
ath hw hh 


WAAAA BKQAAA 
NeEwun wah aAN 


8 
4 
2 
2 
6 
4 
4 
2 
4 
2 


wh pad 


_ 
Onann 
—ww 
mw ND 


40 small chunks 
of concrete 
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CONCRETE FRAGMENT WE:GHTS ANB DIME‘NS'ONS 
FOR 1/24-SCALE SHIELDED REACTOR MONE1L.S 


Model No. 14 


Event Simulated, 160 lb TNT, ! msec, w/Pyrocore 
Vessel 3/4 Full of Water 


Height, Remarks 
(in.) 


Cylinder, 12-in. 


diam. with 
4-1/4-in. diam. 
hole 


Model test Event Simulated, 150 1b TNT, 1 msec, w/Pyrocore 


Vessel 3/4 Full of Water 


Group Weight, | jreight, | Length, | Width, Remarks 
Ib (in.) (in.) (in.) 


Cylinder, with 
4-1/4-in, diam. 
hole at center 


10 small chunks 
of concrete 


SECRET 


Mode} No. 
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CONCRETE FRAGMENT WEIGHTS AND DiMENS:ONS 


16 


Event Simulated, 


210 Ib NT, 
Vessel Full of Water 


FOR 1/24-SCALE SHIELDED REACTOR MODELS 


Tinsec, w/MDF 


Piece es Height. 
Number ae Cin) 
4 
4 
4 
4 
4 
4 
4 
4 
2-1/2 
3 
2 3-3/4 
2 3 
2 3 
ral 2 2-1/4 
2 Z 2-1/2 
22 2 17 1 2-1/4 
12 1412 2-1/2 
13 1] ie 3 
3h ] 9 2-1/2 
14 1 3 2-1/2 
30 1 3 2 
15 byo2 1-1/2 
23 1 
26 te io] 9 2 
32 145 1-1/2 
9 15 1-3/4 
33 12 1-1/2 
27 iW ! 
38 Mi 2-1/2 
24 9 2 
20 8 1-1/2 
28 8 2 
28 B 2 
36 6 2 
37 6 2 
Py 5 2 
43 2-3/4 1-1/2 
46 2-3/4 1-1/2 
44 2-1/4 1 
45 2 1-1/2 
52 1-3/4 2-1/2 
47 1-1/2 1-1/4 
53 | 1-1/2 1-1/2 
418 b-1s4 1-1/2 
419 Il i/2 
40 4 1/2 
51 ya 1/2 
44 v2} By Ss ) 
_ 31 3 
i 
~ \ 6 | 
f 7 | 
| ' 
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| Length, Width, 
Cin.) (in.) 
7 | 10 
1 6 
8 8 
5-1/2 1 9 
8 7 
6 6-1/2 
7 5 
6 4 
5 4-1/4 
4 4 
6 3-1/2 
6 3-3/4 
5-1/2 4-1/4 
5 4-1/2 
5 3-1/4 
5 4-1/2 
4-1/2 3 
3-1/2 3-1/2 
5 2-1/2 
4 2-1/2 
4 3-3/4 
5 2-1/2 
4 2-1/2 
4-1/2 2 
4-1/2 3 
5 2 
3 3 
3-1/2 2-1/2 
3 1-3/4 
3. 2 
3-3/4 2 
2 2 
3 3 
3 2-1/2 
3 2-1/2 
2-1/2 1-3/4 
1-3/4 1-1/4 
2-1/2 1-1/2 
2 i | 
1-3/4 | 1 | 
1-1/4 1 
1-1/4 1-1/2 
1-1/2 } 
1-1/2 3/4 
1-1/2 ae 
2 1/2 
2 ! 
1 : 3/4 
= 
| 
| | 
' | 
| 
hae ete he peas Ty. - Ta sa 


+ 


Dist. of Fr 
from Original Model 


apments 


Position, If > 10 ft 


ft in, 


Remarks 


Approximately 
200 sinall 
chunks of con- 
crete (6 oz. to 
1/2 ov range) 


Approximately 
125 very small 
chunks of cone 
crete (below 

V2 021m sine) 
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CONCRETE FRAGMENT WEIGHTS AND DIMENSIONS 
Mictchnes te FOR 1/24-SCALE SHIELDED REACTOR MODELS 
ode Oe 


Event Simulated, 160 lb TNT, 1 msec, w/Pyrocore 
Vessel Full of Water 


Individual 


Grou 


"Weight, Weight, | eight, | Length, | Width, Remarke 


Pp 


Ib 
145{ 145 Cylinder, 12-in. 

diam. with 
4-1/4-in, diam. 
hole 

4 12 

1 4 1 

1l 1 

7 

4 

2 

1-1/2 

1 15-1/2 
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| Remarka 


Position, if > 10 ft | 


fram Original Model 
26 


Cin 


INT, } msec, w/MDF 
Width, 


1 
4 


Length, * 
Cin) 


| 
L 
: 
| 
| 
! 


Vessel 3/4 Full of Water 


Herght, 


Event Simulated, 210 1b 
( in.) 


19 


CONCRETE FRAGMENT WEIGHTS AND DIMENSIONS 
FOR 1/24-SCALE SHIELDED REACTOR MODELS 


Mode] No. 


chant a af cone 
crete (how to 


Approndravtely 
1/2 07 range) 


TEU worst 


2 
| pea/2 


1-1/2 


wll 
' 


Apprasiia tlety 
20 very 5 


a 


‘ 


{ 
' 
i 


BS 3s Seraene eta tajccai re 


F-7 


SECRET 


APPENDIX G 
CONCRETE FRAGMENT WEIGHTS AND DIMENSIONS 
FOR 1/12-SCALE SHIELDED REACTOR MODELS 


CONCRETE FRAGMENT WEIGHTS AND DIMENSIONS 
FOR 1/12-SCALE SHIELDED REACTOR MODELS 


Energy Source: Pyrocore 
Period: l msec 
Water in Vessel: 100% 


Model No. 
Event Simulated 
TNT (1b) 


Megawatt-Seconds 


Model Wt. (1b) 


% Fragments 

<1lb 

1-51b 1.4 

5 - 10 lb , 1.3 

10 - 15 1b 1.9 

15 - 20 lb 1.4 

20 - 25 lb 8.7 
> 25 lb 83.9 


No. Fragments 
<11b 31(1616) 
1-51b 10 
5 - 10 lb 
10 - 15 lb 
15 - 20 lb 
20 - 25 lb 
> 25 lb 


No. Fragments 
Excluding Fines 


* ec ate : 
Figures in parentheses indicate "Fines". 


SECRET 


0.5 
1.7 
1.1 
3.9 
13.2 
78.6 


22(1050) 
4 
3 
1 
3 
8 
14 


55 


SECRET 


1.4 
3.2 
6.1 
3.6 
8.0 
76.1 


47(99) 
8 


6 
7 
3 
5 


20 


96 
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APPENDIX H 
TRAJECTORY OF AN AIR PARTICLE 
DURING THE POSITIVE PHASE OF A NUCLEAR BLAST 
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APPENDIX H 
TRAJECTORY OF AN AIR PARTICLE DURING THE POSITIVE PHASE 
OF A NUCLEAR BLAST 


Using the 1957 ''The Effects of Nuclear Weapons'' {Ref. 20), the 
velocity-time history of an air particle may be found at any fixed point on 
the ground surface for a surface burst of any yield. If we know the velocity 
of an air particle u, at one location, we can predict its new lccation after 


a short time interval At) (e.g., old location + u AT,). Repeating this 


procedure enables us to establish the horizontal Ee ae relationship 
of the particle. Because so many factors enter into the determination of 
the particle velocity as presented in reference 20, an otherwise straight- 
forward integration becomes a rather involved bookkeeping problem. The 


following outline indicates the steps involved in the numerical integration: 


I, Select 
A. Weapon Yield 


B. Overpressure at front of forest, Py 


Il. Compute (Time: T = 0) 
A. Shock Velocity U; for P) (see Fig. 3.80, Ref. 20) 
B. Particle Velocity u, for p) (see Fig. 3.80, Ref. 20) 
C. Distance from ground zero (GZ) to front of forest, X 
1. Fig. 3.94a; 1-KT surface burst; Ref. 20 
2. Scaling Law; Eq. 3. 86.1 Ref. 20 


Ill. Select T = AT, 


IV. Compute 


A. Distance traveled by particle in AT) Ax, = uy) AT, 


Ay, = U,AT 


B. Distance traveled by shock front in T )47T; 


1 
C. Time gap between shock front and particle: 


t) = u, - u,) AT, /U, 


SECRET 


a 
w 


SECRET 


f 


D. Overpressure at distance x, + Ax, from GZ: PS 


1 
1. Fig. 3.94a; 1-KT surface burst; Ref. 20 
2. Scaling Law; Eq. 3.86.1 Ref. 20 

E. Overpressure at distance x, + Ay, from GZ: P> 


1. Fig. 3. 94a; 1-KT Surface burst; Ref. 20 
2. Scaling Law; Eq. 3. 86.1 Ref. 20 
F. Shock Velocity U; for P> (see Fig. 3. 80, Ref. 20) 


G. Duration of positive phase at x, + Ax): tia 


1, Fig. 3.96; surface burst; Ref. 20 
2. Scaling from Art. 3. 88, Ref. 20 
H. Overpressure behind shock front p(t,) 


1. Compute t)/tyy 
2. Compute p(t)) from Eq. 3. 82.1 (Ref. 20) using P5 


I. Particle Velocity u, for p(t,)5 (Fig. 3.80 Ref. 20) 


V. Select (T = AT,) 


VI. Compute 


Ax, = u,AT, 
Ay, = U,AT, 
t, = [dy, + Ay, - (Ax, + Ax,)|/U, 


P3 atx, + Ax, + Ax, from GZ 


P3 atx, + Ay) + Ay> from GZ 


Bh > Oh 


The curves in the figures referred to in the outline were all fit 
with analytical expressions and the entire procedure was programmed for 
the UNIVAC 1105, Figure H-1 shows the detailed flow diagram of the 


program. 
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d= afw't 3 


} * No. 
LAD oe toes AED 


_o 
Yes Yes 
(328) te (1s i= 8? TCHL)G)I8) 


(Gi) p= Ad B @ 
oO a 
Read Input 


Cs a 
EID ee eae) CRETE = CREED ee REED 
i 4. 
No - x 
1f3 
4 


Yes G:) 
Ca dow dg 


pam }—fow'? Hye 3 
©) No 


In a. 1,0? 


res 


T(33)KM 24) 


pap (1.0- aye" Be (1g +6p/7»,) fe 


T(35)(a)(4. 0) 


oer ¥ : 
uc ape UeG 
AO py 
5 


THON) 
TAMA GMP SM) 
TO Mcayin 
rts) 


ote 
Oo Ponete . Mattiptie tian ao 


Figure H-1 Detailed Flow Diagram For The Particle Velocity Integration 
Se hy ee ration. 
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THE VULNERABILITY OF 


ANTENNA SYSTEMS 
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APPENDIX J 
THE VULNERABILITY OF ANTENNA SYSTEMS 


The severity of the debris hazard to antenna systems 1s only 
one phase of debris problems associated with hardened sites, but it is 
used to describe the estimating procedures that can be applied to debris 


problems in general, using data collected in this report. 


The approach used in this study made estimates in several 
manners and noted their consistency. First, estimates of the maximum 
range of antenna vulnerability were based on the correlation of maximum 
debris distance and equivalent yield observed in Chapter Two. Debris 
distribution data from the Pantex Ordnance Plant event were then scaled 
up to nuclear yields. This was done by scaling up the ground ranges 
according to the cube-root-of-yield scaling, although it may have been 
more appropriate to also scale up the total volume of fragments in a manner 
which accounted for the total material volume of the nuclear crater. The 
result was perhaps optimistic, resulting in fewer fragments. Next, debris 
environment was estimated independently by scaling up DANNY BOY findings 
to high yields. This was done by scaling both ground ranges and fragment 
densities according to cube-root-of-yield scaling. A third estimate of 
debris environment was based on use of the hydrodynamic model of crater 
formation developed by Brode and Bjork at RAND Corporation’. An analytical 
solution to the debris environment was obtained by considering peak 
velocities from the hydrodynamic model to be initial fragment velocities, 
and following the trajectories of the ejecta from the crater to ultimate 
impact with the ground. Results of these three approaches are consistent 


in predicting severe debris environments under likely hardening criteria. 


"H. L. Brode and R. L. Bjork, ''Cratering from a Megaton Surface 
Burst", Paper L, Proceedings of the Geophysical Laboratory - 
Lawrence Radiation Laboratory Cratering Symposium, UCRL-6438, 
University of California, October 1961. 
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There are two potential sources of debris hazards to antenna 
systems: throwout material from the crater and loose material or broken 
structure picked up and transported by the blast winds. The second source 
of debris can be eliminated by clearing an area around the antenna. No 
simple defense against crater throwout debris is apparent. This problem 
is therefore restricted to prediction of the debris hazard resulting from the 


crater throwout alone. 


A comprehensive method of analysis of debris effects emanating 
from nuclear detonations has yet to be developed. To a great extent the 
character of the hardened sites is responsible for this state of affairs. The 
hardened missile silo with a reinforced concrete cover as the only exposed 
component was considered relatively invulnerable to debris damage. The 
sole problem was to predict the weight of debris on the closure after an 
attack so that sufficient power could be designed into the closure operating 
mechanism. Communications systems, however, pose entirely different 
problems. Components of antenna systems, for example, may be quite 
vulnerable to debris damage because of their electrical essentials. Thus, 
the hardened antenna design must be based on such criteria as density, 


size, and energy of debris particles to be expected at the antenna location. 


Some studies which border on,or are corollary to,the problem 
have been performed. These include studies of missiles from accidental 
explosions (Chapter Two), studies of debris distribution from nuclear tests 
(JANGLE U and DANNY BOY), and analytical studies dealing with the crater 
formation problem (RAND). Each of these is used to estimate the debris 


hazard for the hypothetical antenna structure. 


It should be emphasized that at present, it is only possible to 
assess the debris problem very roughly. Our intention is to objectively 
review the available data and to compare predictions based on different 
source material. Hopefully, the results will be consistent. We do not 
contend that this study completely settles the debris question for hardened 
sites. We do feel, however, that it is the most reasonable approach to 
the problem within the current state-of-the-art, -- i.e., short of further 


nuclear testing. 
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J-1 Estimates Based on High-Explosive Debris Data 


The Armed Services Explosive Safety Board data on 206 HE 
detonations ranging in magnitude from 8-1b of tetrytol to 9, 000, 000 lb of 
ammonium nitrate were collected (Chapter Two). A statistical analysis 
relating the maximum missile distance to the weight of explosive (TNT 
equivalent) was made. On a log-log plot a linear regression line to best 


fit the data was found to be of the form, 


log, pPay = 2.950 + 0.322 log, W (J-1) 
where 
Dy = maximum missile distance, ft 
W = equivalent weight of TNT, lb. 
Similarly, a quadratic regression line to best fit the data was found to be: 
2 
10g) Day = 2.960 + 0. 347 log) pW - 0.016 (log) yW) (J- 2), 


Note that Eq. (J-1) indicates one should scale according to w°: aoe A 


standard error was found to be 2.47 Dy: Equations (J-1) and (J-2) were 
applied to the antenna vulnerability problem using an equivalence factor of 
0.50 to relate nuclear yield to TNT equivalent; the results are presented in 


Table J-1. 


Table J-1 
LIMIT OF ANTENNA VULNERABILITY BASED ON HE DATA 


Based on Based on 
Weapon Linear Regression Line Quadratic Regression Line 
Yield, Maximum Range of Maximum Maximum Range of Maximum 
(MT) Debris Debris Distance for Debris Debris Distance for 
Distance, One Standard Error Distance. One Standard Error 
(miles) (miles) (miles) (miles) 


7.9to 48.0 6.4 2.6 to 15.7 
9.8 to 59.9 


20 30.3 12.3 to 75.0 7.6 
50 40.7 16.5 to 100.7 8.4 
100 50.9 20.6 to 121.0 9.1 
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It is interesting to note that a data point corresponding to 5, 200, 000 lb of 
TNT (2.60-KT nuclear weapon) was included and that the maximum missile 
distance for this point was 3.5 miles. Thus, we extrapolated four cycles 


from six cycles of data. 


The relationship between the high explosion debris problem and 
the debris emanating from the crater of a nuclear weapon is certainly ques- 
tionable. The HE debris generally results from buildings and equipment 
in the immediate vicinity of the explosions. The mechanism by which this 
debris is formed is different from mechanisms by which a crater is formed. 
Nuclear detonation is accompanied by rather substantial winds whereas 
in the HE explosion these winds are essentially absent. Nevertheless, it is 
desirable to take into consideration the large body of data that is available 
for the HE debris problem, particularly in view of the scarcity of nuclear 
test data. As a matter of fact, since essentially all debris associated with 
an HE detonation originates near the point of detonation, one can argue that 


this debris is, in fact, similar in origin to throwout debris. 


The second item of interest is the distribution of debris outward 
from the point of the explosion. A detailed debris study which included 
complete descriptions and final locations of the debris was carried out at 
Pantex Ordnance Plant (Ref. 11). The explosive was 2000 lb of TNT 
detonated in a reinforced concrete bunker. For application to the antenna 
vulnerability problem, all reported ranges were scaled up by the cube-root- 
of-yield law. The total volume of fragments at the Pantex study was 
estimated to be 1000 cu ft, whereas the crater volume for, say a 20-MT 
weapon surface burst is about 3.7 x 10? cu ft. Rather than using the ratio 
of 3.7 x 107/107 to scale the number of fragments, we used the more 
optimistic (resulting in fewer fragments) cube-root-of-yield factor (215 
for this case), The resulting missile density as a function of ground range 
is shown on Fig. J-1. The dashed portions of the curve are extrapolated 
to expected antenna locations. These results make it apparent that. debris 


problems are critical for antenna systems. 
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J-2 Estimates Based on Nuclear Test Results 


Armour Research Foundation (Ref. 2) recently completed an 
experimental study of crater throwout from an underground nuclear detona - 
tion; these results are derived completely from that work. This study was 
part of Project 7 DANNY BOY, a 0.43-KT nuclear device buried at a depth 
of 110 ft in basalt on the Buckboard Mesa at the Nevada rest site. The ex- 
pected crater zone was salted with a variety of objects which were located 
after the blast. In addition, the natural debris was tabulated in certain 
areas. Of particular interest here are the natural debris radial distribution 
charts reproduced from reference 2 in Fig. J-2, J-3, and J-4. AreaslI, I, 
and III in these figures correspond to different orientations with respect to 
ground zero. The results must be scaled up to the cases of interest here. 
The ground ranges are reasonably scaled by the cube-root-of-yield law. 
Scaling of the fragment density expressed as fragments~per-square-foot 
requires some discussion. If the total number of fragments is assumed 
proportional to the crater volume, then the number of fragments scale directly 
as the weapon yield (crater diameter and depth each scale as the cube root 
of yield). The area over which these fragments are distributed scales as 
length square or as yield to the two-thirds power. Therefore, the fragment 
density scales as w/w2/3 or as cube-root of yield. The results of 
Fig. J-2, J-3, and J-4 for the three sectors are averaged and scaled 
up to a 1-KT weapon as a standard. The result is shown on Fig. J." 

This is then used to determine debris density as a function of range for any 


weapon yield in kilotons. 


This has been done for an antenna system location of 4.0 miles 


for selected weapon yields. Results are shown in Table J 2. 
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Figure J-2 Throwout Density, Area I, DANNY BOY Event. 
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Figure J-3 Throwout Density, Area II], DANNY BOY Event 
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Table J-2 
DEBRIS DENSITY FOR ANTENNA SYSTEM LOCATION 
AT 4.0 MILES GROUND RANGE 
(Based on extrapolation of DANNY BOY data) 


Debris Density, (fragments/sq ft) 


l in. to 6 in. 6 in. to 12 in. 12 in. to 18 in. 
Equivalent Diameter Equivalent Diameter Equivalent Diameter 


(MT) (0. 0003 to 0. 065 ft?) (0.065 to 0. 523 ft°) (0. 523 to 1.77 ft”) 


0. 86 0. 00037 

43.0 0.171 
1, 260. 0 41.3 
11,400. 0 1,410.0 


Based on the assumptions made with regard to scaling, it is obvious from 
Table J-2 that debris would pose a major problem. Note that particles of 
the size shown in Table J-2 were found in the DANNY BOY test. Thus, 

in applying the results of Table J-2 to antenna vulnerability problems, all 
three sizes of particles must be considered simultaneously. No attempt 
was made here to scale" particle size which probably varies with both 


weapon yield and soil type. 


These results are more severe than those predicted from HE 
data (Fig. J-2). Recall, however, that the number of fragments was 
scaled up by only the optimistic cube-root-of-yield scaling. Therefore, 


we feel the Table J-2 results are more significant. - 
J-3 Estimates Based on Analytical Studies 


A completely analytical treatment of this problem depends 
primarily on a mathematical model to treat the crater formation problem. 


Such an analysis gives initial velocity vectors of material or particles leav- a 


ing the crater, and the subsequent motion of these particles can be followed 
by means of standard trajectory analysis. We use a crater model devised Se 
at the RAND Corporation and then compute debris density at ranges from 


ground zero. The RAND model is first discussed, the trajectory analysis 
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is then presented, and finally, results are presented which are applicable 


to antenna systems. 
J-3.1 Application of the RAND Crater Model 


Brode and Bjork studied the formation of.a crater resulting 
from a 2-MT weapon surface burst, assuming the material in the crater 
zone to be rock (tuff). Their interest was primarily in the early period of 
the crater formation when the pressure acting on the rock medium is very 
much greater than the shear strength of the rock. They, therefore, 
assumed the hydrodynamic model valid and numerically integrated the 
appropriate field equations. Pressure and velocity fields in the crater zone 
are presented. The velocity fields are reproduced here as Fig. J-6 
through J-9. 


The RAND model is formulated in Eulerian coordinates so 
that the velocity vectors represent the velocity of the mass currently at the 
point in space indicated by the base of the vector. It is therefore not 
possible from the available data to rigorously follow the motion of a 
specific mass of crater material. Rather we used these data by assuming 
that the peak velocity at each point in the crater is the initial velocity of the 


mass at that point. 


A grid was established to roughly cover the crater and the data 
of Fig. J-6 through J-9 were used to determine peak velocities at each 
point of grid. The grid is shown in Fig. J-10 with grid points identified. 
The velocity vectors which are most severe from the viewpoint of throwout 


were selected; the resulting velocities are shown in Table J- 3. 
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Figure J-6 Velocity Field in Crater (T= 0.1026 msec) 
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Table J-3 
INITIAL TABLE VELOCITY 


Grid Point (see Initial Horizontal Initial Vertical 
Fig. J-10 for Lo- Velocity, Velocity, 
cation in Crater) (meters /msec) (meters /msec) 


OWDND UhPWN eH 


2 99929 SSee9 
2 99°99 Se9e° 


The velocity data in Table J-3 were used as the initial velocities 
in computing debris particle trajectories. Because the hydrodynamic model 
does not predict debris fragment size, trajectories were computed for a 


range of particle sizes. 
J-3.2 Trajectory Analysis 


Consider the motion of a particle through a medium such that the 
drag force acting on the particle 1s proportional to the square of the rela- 
tive velocity between the particle and the air. It,is assumed that the verti- 
cal and horizontal motion of the debris particle are decoupled. This is true 
if the center of pressure of the particle coincides with its centroid for all 


orientations so that no rotation occurs. 


The equations of motion are then 


s= KOP (u- xy? (J- 3) 
di! GRU, 2? 
yr Sy tg wo) 
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where 
x = horizontal coordinate 
y = vertical coordinate, positive downward 
('; = differentiation with respect to time 
a = aerodynamic coefficient, (projected area x drag coefficient) —-~ drag coefficient) 
p = mass density of air 
u = air particle velocity 
g = gravitational constant 
ie -{* 1; u>x 
-1; ucx 
ig t liys<o 
-l,y>o 


These are Riccati-type nonlinear differential equations, and can be linear- 


ized by a simple transformation of coordinates. 


The horizontal equation of motion (J-3) can be linearized by 


the substitution, 


a 2. 58 
“ " kap s- (J -5) 
Note that x =- E lo s 
~ kap Be * 
thus, 
oe . Ke a” 2 2 
Ss + kapus + ———f u’ sS= 0. (J -6) 


4 


Equation (J-6) is a linear differential equation with variable coefficients 
because of the variation of particle velocity u and air density p with time. 
This could be numerically integrated but, because of our interest in 
relatively large times, the computation time would be prohibitive, and the 
resulting cumulative error undoubtedly sizable. Equation (J-6) is there- 
fore solved by assuming that u and p are constant over an interval of 


time. The solution can then be extended in time by matching initial 
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conditions after each interval of time, and then changing u and p toa 
new constant value for the next interval. If u is constant, Eq. (J-6) 
has a solution 
a ~ Yapk ,. 7 
s = (C, +C, t) exp]- —j- t (J -7) 
where 


C)> and Cc, are constants of integration. 


The vertical equation-of-motion Eq. J-4 can be linearized by 


the substitution, 


2 s 
YS ep ee ee) 
Then | 
2 : 
yes gdp longs eG 


Then reduced equation then becomes, 


Z +k! “88. zZ=0 (J-10) 


Recalling that k'= +1 for ygo and k' = -1 for y>o, Eq. J-10 has the 


solution, 


z= C , cos\/ 88 t+C, sin\¥ (J=11) 


for v< o (i.e., on way up) 


z= c exe [Vee | t+ op | V5 | 


for y>o (i.,e., on way down). 


Equations (J-11) are used to compute the total time of flight for the particles, 
and then the total horizontal distance traveled is determined from the value 


of s, as computed from Eq. (J-7) at the time the particle hits the ground. 
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J-3.3 Numerical Results 


The simultaneous solution of Eq. (J-6) and (J-11) was carried 
out on the UNIVAC 1105 digital computer. The initial velocities presented 
in Table J-3 were used as initial conditions and solutions were obtained for 


l-in., 6-in., and 12-in. diameter particles. 


The analytic forms of the weapon parameters used in the 
computer program are taken from reference 20. The actual values as 
taken from the computer program for a 20-MT weapon are plotted in Fig. 
J-11. For lack of better data, the close-in values for overpressure shock 
velocity and particle velocity were taken to be constant from ground zero out 


to a scaled ground range of 100 ft. 


Solutions were obtained for particle velocity and air density 
assumed to be constant in 0. 25-sec and 0. 1-sec intervals. The numerical 
results differed by less than 2 percent so the computer runs were finally 
made using the 0. 25-sec interval. Numerical results for flight time, hori- 


zontal distance traveled and final velocity are given in Table J -4. 


These data were then converted to fragment density values. 


Consider the crater to be broken up into annular rings as shown in Fig. J-10. 


The material in each of the three horizontal layers was first distributed 
over the impact area. This was done by assuming that the material from 
each ring is spread at constant depth over a radial distance equal to the 
difference in the computed maximum trajectory distance of adjacent points. 


The fragment density in that region is then given as 


ip eae —y—+ (3-12) 
2 T +“ mar_-r 
3 f ° i 

where 

N = fragment density 

V = volume in crater as given in Fig. J-10 

T,. = fragment radius 

rT, = outer radial distance for ring of interest 

r. = inner radial distance for ring of interest. 


SECRET 


J-22 


‘uodeam LW-0z7 40} tojourereg uodeam TI-¢ eandty 


(43) ‘eduey punory 


SECRET 


000‘g1 000‘91 000‘FI O00'ZI 000 ‘OI 0008 0009 000% 0002 0 
Ooo! 


exereys 

| 'd 
o} 
a 
000S ey o 
< 

- : 

! U QQ 

eae { Z ! 

n i KR 
io) | “ 
S) OO0L pase e 
| 
= 5 
8 | E. 
0006 }---——--5 # 
Ss | ae 
it) 
o 
Ass 

Soar 
000 ‘II -— -~ -+ -~~-- 4306 see 
2? 
—_— 
| =~ 
rr) 


Sean Ee ss? 2 S| | aemateal nal’ 


oe oe 


~~ 
‘ 
a 
(298) (sd3)  (s@qtur) (sdz) (seqrur) | (Sas) (sdy) = (Saqyur) | (uotzeD07 103 
‘eauity, ‘AoOoTaA ‘uOTZBDOT ‘AQIDOTaA ‘UOTVeD0T | fourty ‘AOOTaA ‘uOT}yeDOT | OI-£ *3taq 2a9) 
ws Teurd Teura [94stta yeurdy reura {3st {eur st Teut a taIeID 
SaTI1W1eq YOuy~aaTany seToaeg your-xtg S9TIIWAeq YOuT-2uo Uy wrod — 
NOdVGM LW-02 YOU VLVd YOd SISATVNV AYOLOALVUL — 


p-f eTdeL 


ke, 


freee 4b 


SECRET 


The effects of the three layers were then added. 


The results of Table J-4 were then transposed to fragment. 
density-distance curves and, for consistency, scaled down to a 1-KT weapon 
(both the range and fragment density were scaled as the cube root of yield). 
The result is plotted on Fig. J-12. It should be emphasized that the particle 
sizes on Fig. J-12 are not predicted by the model. Figure J-12 represents 
three solutions assuming that all of the crater material breaks up into the 
Same size particles. Therefore, only one of the three curves in Fig. J-12 
should be used at one time. This can be compared with the DANNY BOY 
results in Fig. J-5. The dashed line drawn in Fig. J- 12 is from the 
DANNY BOY results, assuming that all debris fragments break up into 
one-inch radius particles. It is not surprising that the RAND model gives 
more severe results because of the assumed hydrodynamic behavior of 


the soil and the neglect. of in-flight collisions between debris particles. 


An interesting result was obtained while studying the trajectories 
of the debris leaving the crater. The total flight distance was found to be 
independent of the initial horizontal velocity component for the range of 
horizontal velocities predicted by the RAND model. In other words, the 
horizontal motion 1s determined completely by the blast winds for large 
weapons. The initial vertical velocity, of course, determines the time of 
flight upon which the total throw of a particle is very dependent. Based on 
this brief analysis, 11 appears that for a megaton-yield weapon, the surface 
burst (or perhaps partially buried burst) results in the most severe debris 


problem. 
J-4 General Consistency of Results 


Crater throwout. debris was studied with particular emphasis 
on predicting its severity to antenna systems. The state-of-the-art is such 
that a completely reliable evaluation of this problem was not possible. The 
objective here was to examine available experimental data and analyses ,. 
to establish bounds on the magnitude of the problem for antenna systems. 
First, the voluminous data which exist on debris resulting from high 
explosive detonation -w:e studied. Results are presented in Fig. J-1 and 


Table J-1. It can be seen that the maximum missile distances predicted 
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are beyond likely antenna locations, and that at expected antenna locations 
the debris density is quite high. Secondly, the results of the DANNY BOY 
nuclear test are extrapolated to make predictions for antenna systems. It 
is assumed that both distance and fragment density scale as the cube root 
of yield. The resulting fragment density-distance curve for a 1-KT weapon 
is shown in Fig. J-5. This can be readily scaled up to real antenna 
situations and againa severe debris problem is predicted. Third, a 
completely analytical solution was obtained by following the motion of a 
particle from the crater to impact with the ground. Brode's hydrodynamic 
model was used for the crater formation phase of the motion. The resulting 
debris density-distance relationships, which indicated a relatively severe 


debris problem for expected antenna locations, are shown in Fig. J-12. 


Therefore, while no conclusive evaluation of the problem was 
possible, all available data point to a very critical debris problem for 
antenna systems. There are two possible approaches to the design of 


antenna systems insofar as the debris problem is concerned. 


First, locations could be restricted to those where the debris 
density was at a specified low level. Secondly, some degree of hardness 
could be provided against debris particles. ‘At first look this would seem 


a most challenging task. The debris will have terminal velocities very 


close to the particle velocity, which can well be over 1000 fps. At these 


velocities it seems apparent that even quite small debris particles would 


be capable of damaging any fragile critical elements. 


Certainly much has yet to be learned about the crater throwout 
problem. The soil type, depth of burial, and yield must all have some 
effect. None of this is understood at present. The results presented here 
are all based on a rock-like crater material. The debris particle size 
must certainly be a function of soil type. A sandy material should produce 
very small fragments that would tend to sand blast rather than fracture 
antenna elements. Also, as can be seen from Table J-4 sand would not be 


transported as far as a soil which breaks into large fragments. 
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The extremes of depth of burial were included in this study. 
DANNY BOY represents a completely buried shot at optimum depth* whereas 
the Brode crater model is for a surface burst. As remarked earlier, it 
appears that the surface burst will give the more critical debris problem 


for megaton range weapons. 


The dependence on weapon yield (scaling) represents perhaps 
the most important unknown. Most of our results depend on an assumed 
cube-root-of-yield scaling law, although the analytical results were only 


scaled up from a two-megaton weapon. 


Perhaps the most surprising result of this study is the 
consistency of the results of the three methods studied when so many unknown 
factors exist. This leads to some degree of confidence in the predictions 


of the study. 


* Optimum depth of burial is that depth which produces maximum volume 
of apparent crater. It results in ejecta having trajectories with 
pronounced vertical components. 
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